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PREFACE 


Choosing  a  beachhead  site  for  HF  D/F  in  support  of  an  amphibious 
operation  usually  involves  a  choice  among  mediocre  or  even  poor  loca¬ 
tions.  The  number  of  suitable  sites  is  restricted  by  both  the  need  for 
security  and  the  limitations  imposed  by  D/F  requirements  themselves. 
The  latter  requirements  may  be  broken  up  into  two  general  categories: 

a)  Physical  characteristics  (i.  e.  ,  is  it  large  enough  and 
sufficiently  clear  and  level  for  antenna  erection?  )  and 

b)  Electrical  characteristics  (i.e.  ,  are  there  nearby 
scattering  objects  such  as  hills,  equipment,  etc.?). 

As  severe  as  these  limitations  are,  there  is  at  present  no  formal  site 
classification  scheme,  nor  a  tactically  usable  method  of  electrical  site 
evaluation.  It  is  clear  that  the  electrical  characteristics  must  be 
determined  before  any  meaningful  classification  can  be  made. 

The  purpose  of  this  report  is  to  present  a  quantitative  basis  for  site 
classification  (see  Table  P-I)  and  a  significantly  simplified  method  of 
measuring  the  electrical  characteristics  of  certain  sites.  The  body  of 
the  report  deals  with  the  method  of  measurement.  This  method  seems 
highly  promising  on  a  theoretical  and  scale -model  experimental  basis, 
but  for  final  validity  awaits  spot-checking  at  HF  frequencies. 

Table  P-I  gives  some  examples  of  the  physical  sizes  required  for 
some  given  antennas.  For  example,  the  Wullenweber  requires  a 
nearly  1,000-foot  diameter  and  the  AN/TRD-21  requires  75  feet  just 


PHYSICAL  PARAMETERS 


Clear  Diameters  (in  feet)  Examples 


300  -  1000 

Wullenweber 

100  -  300 

Doppler 

30  -  100 

Adcock 

10  -  30 

Elevated  H 

3  -  10 

Whip 

ELECTRICAL  PARAMETERS 

si 

Examples 

3-10"3  -  10'3 

Null  depth  of  28  - 

33  dB 

co 

i 

o 

i 

co 

i 

rsj 

i 

o 

Null  depth  of  23  - 

28  dB 

-2  -2 

3-10  -  10 

Null  depth  of  18  - 

23  dB 

CNJ 

1 

o 

1 

CO 

1 

1 

o 

Null  depth  of  1 3  - 

18  dB 

OJ 

o 

1 

*— • 

o 

1 

Null  depth  of  8  - 

13  dB 

TABLE  P-I  PHYSICAL  AND  ELECTRICAL  PARAMETERS 


for  the  antenna  elements.  Examples  are  also  given  for  the  electrical 
2 

parameter,  k  ,  which  is  a  ratio  of  average  scattered  power  to  direct 

power  This  parameter  was  related  to  null  depth  of  a  figure-eight 

pattern,  not  only  to  make  the  classification  less  arbitrary,  but  also  to 

give  the  user  more  of  a  "feel"  for  i*r  meaning.  A  null  depth  of  20  dB 

may  be  considered  average  for  most  siting  conditions,  while  30  dB 

indicates  a  good  location  and  a  minimum  of  scattering  obstacles. 

Similarly,  a  10  dB  null  depth  would  indicate  a  relatively  poor  D/F 

2 

site.  The  method  of  converting  k  into  equivalent  null  depth  is  given 
in  Appendix  A. 

The  currently  available  method  of  determining  the  electrical  character 
of  a  site  is  extremely  complicated  and  time-consuming  under  the  best 
of  circumstances.  In  order  to  test  a  D/F  site,  one  must  move  a  trans¬ 
mitter  around  the  receiving  antenna  at  a  given  distance  and  take  readings 
for  each  azimuthal  position  at  all  frequencies  of  interest  The  disad¬ 
vantages  here  are  obvious.  Keeping  a  known  distance  from  the  antenna 
of  several  thousand  yards  is  impossible  under  beachhead  constraints. 
Even  if  it  were  feasible  to  maintain  a  fixed  distance,  the  allowed  azi¬ 
muthal  positions  would  fall  far  short  of  those  required. 

The  new  technique  described  in  this  report  will  enable  a  transmitter 
at  one  or  two  fixed  positions  to  scan  the  frequencies  of  interest  and 
extract  the  same  information  as  is  obtained  by  the  more  cumbersome 
procedure  Chapter  4  of  this  report  gives  the  general  experimental 
techniques.  A  practical  guide  to  the  field  use  of  this  method  may  be 
summarized  as  follows: 


xiii 


1)  Place  a  target  transmitter  approximately  100  wave¬ 
lengths  distant  from  the  site. 

2)  Sweep  the  transmitter  over  a  one -octave  band 
centered  on  the  frequency  of  interest. 

3)  Continuously  record  the  received  field  strength 
(measured  with  a  whip)  at  the  site  as  a  function 
of  frequency. 

4)  Construct  a  moving  average  field  strength  by 
averaging,  at  each  frequency,  the  field  strength 
over  a  0.  1  octave  band  centered  on  that  frequency. 

5)  Take  the  ratio  of  the  measured  field  strength  to 
the  moving  average  field  strength  at  40  equispaced 
points  over  the  octave  band. 

6)  Determine  the  variance  of  these  40  samples.  This 

2 

is  the  value  of  k  . 

The  actual  numbers  given  are  somewhat  arbitrary  and  other  choices 
may  be  found  to  be  satisfactory.  Steps  4  through  6  may  be  quickly 
performed  with  a  special  portable  device.  While  this  device  is  not 
yet  available,  a  diagram  is  shown  in  Figure  P-1. 

Having  obtained  the  electrical  parameter,  it  is  feasible  to  assign  a 
physical  parameter  on  the  basis  of  available  site  diameter  as  shown 
in  Table  P-I.  For  conciseness,  we  may  assign  coded  labels  to  all  of 
the  parameters  in  Table  P-I  and  regroup  them  as  shown  in  Table  P-II. 
On  the  basis  of  this  table,  one  may  now  speak  of  a  D/F  site  as  belong¬ 
ing  to  one  of  25  categories ;  e.  g.  ,  A-3,  B-l,  D-5,  etc.  This  removes 
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FIGURE  P-1  A  POSSIBLE  AUTOMATIC  METHOD  OF  ELECTRICAL 

CHARACTERIZATION  OF  HF  D/F  SITES 
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.1 

.1 
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Physical 

Electrical 

n 

Code 

Diameter  (feet) 

Code 

LI 

0 

A 

300 

-  1000 

10-3 

-  3-10'3 

1 

B 

100 

-  300 

3  - 1 0  ~  3 

-  io'2 

2 

n 

C 

30 

-  100 

io-2 

-  3  -  1 0  2 

3 

D 

10 

-  30 

3-10"2 

• 

-  io"1 

4 

E 

3 

-  10 

10-' 

-  3  -  1 0  ”  1 

5 

TABLE  P-II  SITE  PARAMETERS 


the  use  of  such  qualitative  terms  as  "good"  and  "poor.  "  Although  each 
category  defines  a  range  of  values  for  each  parameter,  they  are  well- 
defined  and  subject  to  experimental  verification. 

The  work  described  in  this  report  covers  the  theoretical  justification 
of  the  model  used  and  provides  experimental  verification  for  a  scaled- 
down  model.  The  materials  used  were  chosen  to  correspond  to  actual 
conditions  at  HF.  Ultimately,  the  scheme  will  require  full-scale 
verification  in  the  HF  region,  as  well  as  the  collection  of  enough  data 
to  further  verify  the  statistical  nature  of  the  study. 

JAMES  CARMICHAEL 


xvu 


A  STATISTICAL  ANALYSIS  OF  SITING  ERROR 


IN  RADIO  DIRECTION-FINDING 


by 

Laurence  Anthony  Lyons 


A.  B.  ,  Harvard  University,  Cambridge,  Massachusetts,  I960 
M.  E.  E.  ,  Syracuse  University,  Syracuse,  New  York,  1964 


Dissertation 


Submitted  in  partial  fulfillment  of  the  requirements 
for  the  degree  of  Doctor  of  Philosophy  in  Electrical 
Engineering  in  the  Graduate  School  of  Syracuse  University, 

May,  1967 


APPROVED: 

DATE: 


1 


ACKNOWLEDGMENTS 


My  interest  in  radio  direction-finding  evolved  into  the  work  presented 
here  through  many  invaluable  discussions  with  my  advisor,  Dr.  A.  T. 
Adams.  I  am  deeply  indebted  to  him  for  his  guidance  and  encourage¬ 
ment. 

I  wish  to  express  my  thanks  to  the  management  of  the  Syracuse 
University  Research  Corporation  and  particularly  to  Dr.  B.E. 
Simmons,  for  the  support  given  my  efforts. 

My  wife  has  been  a  continual  source  of  strength  and  confidence  during 
the  completion  of  this  study. 

I  am  grateful  to  the  secretarial  staff  of  the  Syracuse  University 
Research  Corporation  and  Mary  Jo  Fairbanks  of  the  Electrical 
Engineering  Department  for  their  assistance  in  preparing  this  manu¬ 
script.  Mr.  A.  H.  Berical  of  the  Syracuse  University  Research 
Corporation  provided  much  help  in  the  mechanical  design  of  the 
experimental  apparatus. 

This  work  was  sponsored  by  the  Office  of  Naval  Research  and  the 
U.S.  Marine  Corps  under  Contract  NOnr  4798(00). 


3 


ABSTRACT  OF  THESIS 


Reradiation  from  various  parts  of  the  site  is  a  significant  cause  of 
direction-finding  error  in  the  HF  (1  -  30  MHz)  band  Previous  inves¬ 
tigators  have  considered  the  errors  introduced  by  one  or  two  localized 
reradiators  (e.g.  ,  towers,  small  objects).  Only  a  small  minority  of 
sites  is  characterized  by  the  presence  of  such  obstacles  In  this  work 
a  statistical  model  of  a  site  is  devised  to  analyze  the  effect  of  extended 
reradiators  (e.g.,  hills,  buildings)  or  D/F  accuracy.  Such  a  model 
applies  to  the  majority  of  D/F  sites  encountered  in  practice. 

In  the  model  a  series  of  assumptions  leads  to  an  explicit  expression 
for  the  probability  density  distribution  of  the  electric  field  amplitude 
at  the  direction-finder  as  each  of  two  quantities  is  varied: 

A)  the  azimuthal  position  of  a  target  transmitter; 

B)  the  frequency  of  a  target  transmitter. 

If  the  derived  density  distribution  is  expressed  as  a  function  of  the 
ratio  of  the  total  electric  field  amplitude  (scattered  waves  plus  direct 
wave)  to  the  direct  wave  amplitude,  the  distribution,  in  each  case, 
becomes  a  one  -  par  ameter  Rice  distribution  It  is  suggested  that  this 
parameter  be  employed  as  a  single -number  criterion  for  evaluating 
D/F  sites.  The  parameter  which  is  the  ratio  of  the  mean  square 
power  in  the  scattered  waves  to  the  power  in  the  direct  wave,  may  be 
determined  by  calculating  the  variance  of  a  number  of  samples  of  the 
received  electric  field  amplitude  (normalized  to  the  direct  wave). 

These  samples  may  be  taken  over  angle,  with  the  transmitter  frequen¬ 
cy  fixed,  or,  equivalently  over  frequency,  with  the  transmitter 
position  fixed 
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The  validity  of  the  statistical  model  was  examined  theoretically  and 
experimentally.  In  the  theoretical  work,  the  case  of  scattering  from 
perfectly  conducting  hemispherical  bosses  on  a  conducting  plane  was 
found  to  be  well  described  by  the  statistical  model.  In  the  experimen¬ 
tal  work,  the  case  of  scattering  from  irregularly  shaped  obstacles  on 
an  imperfectly  conducting  earth  was  also  adequately  described  by  the 
model. 

The  new  results  contained  in  this  work  are  as  follows: 

1)  The  application  of  statistical  methods  to  D/F  site  selection. 
Although  the  theory  employed  has  appeared  in  the  literature 
on  ionospheric  propagation  and  rough-surface  scattering,  the 
assumptions  leading  to  the  theory  had  to  be  justified  and  re¬ 
interpreted  in  the  context  of  this  problem. 

2)  The  introduction  of  a  s ingle -parameter  evaluation  of  a  site. 
Current  methods  of  site  evaluation  are  rule-of-thumb  proce¬ 
dures,  which  do  not  possess  a  firm  analytical  basis.  The 
single  parameter  proposed  is  a  physically  acceptable  measure 
of  site  quality  and  can  be  directly  related  to  the  ac'.uracy  of 

a  specific  direction-finder  on  a  given  site. 

3)  The  equivalence  of  measurements  at  various  angles  of  arrival 
(A,  above)  to  measurements  at  various  frequencies  (B,  above). 
This  is  probably  the  most  useful  single  result  of  the  investi¬ 
gation.  The  IEEE  recommends  a  rough  technique  for  site 
selection,  which  depends  on  an  examination  of  the  fluctuations 
in  the  received  power  as  a  target  transmitter  is  moved  in  a 
circular  path  about  the  site.  For  large  radius  circles,  this 
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is  generally  a  difficult  procedure.  The  simplified  procedure 
recommended  in  this  work  involves  examining  the  fluctuations 
in  received  power  as  the  frequency  of  a  fixed  position  target 
transmitter  is  varied. 

4)  A  means  of  correcting  for  the  presence  of  a  few  obstacles  in 
the  vicinity  of  the  direction-finder.  The  means  suggested 
becomes  cumbersome  for  more  than  three  obstacles,  but  it 
is  an  improvement  over  present  methods,  which  can  correct 
for  only  one  obstacle. 
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Chapter  1 


INTRODUCTION 

Radio  direction-finding  has  been  investigated  for  over  seventy-five 
years.  The  earliest  known  direction-finder,  the  rotatable  loop,  was 
first  used  by  Hertz  in  1888.  In  various  forms,  it  is  still  in  use  today. 
Since  the  principle  of  its  operation  and  some  of  its  limitations  are 
common  to  many  modern  D/F  systems,  it  is  perhaps  worthwhile  to 
describe  it  briefly. 

The  small  loop  is  a  directional  antenna.  When  it  is  used  as  a  receiv¬ 
ing  antenna,  it  may  be  employed,  in  theory,  to  determine  the  direction 
of  arrival  of  a  plane  wave  of  arbitrary  polarization  by  orienting  it  so 
that  no  voltage  is  induced  in  it.  The  plane  of  the  loop  then  coincides 
with  an  equiphase  surface  and  the  axis  of  the  loop  is  parallel  to  the 
Poynting  vector.  This  property  of  the  loop  caused  it  to  be  used  at 
first  as  a  direction-finder  for  vertically  polarized  ground  waves.  The 
loop  was  mounted  with  its  axis  parallel  to  the  ground  and  was  rotated 
about  a  vertical  axis  until  a  null  was  achieved.  Although  the  early 
loops  did  not  perform  as  well  as  has  been  described,  because  of  inter¬ 
actions  between  the  antenna  and  the  feed  structure,  later  loops  did, 
after  the  introduction  of  the  ba'1  meed,  shielded,  loop.  * 

A  much  more  serious  limitation  to  the  performance  of  the  loop,  and 
most  D/F  systems,  comes  from  multipath  propagation.  This  effect 
was  first  noticed  when  the  loop  was  sufficiently  far  from  the  trans- 
mittei  to  receive  both  ground-wave  signals  and  sky-wave  signals 
reflected  from  the  ionosphere.  Since  the  phenomenon  was  more  notice - 
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able  at  night,  because  of  the  lowered  height  of  the  ionosphere,  it  be¬ 
came  known  as  "night  effect.  "  It  arises  in  the  following  way:  A  loop 
mounted  with  its  axis  parallel  to  the  ground  will  null  on  a  vertically 
polarized  sky  wave,  when  the  axis  of  the  loop  lies  in  the  plane  of  inci¬ 
dence  (the  vertical  plane  containing  the  wave  normal).  It  will  not  null, 
however,  when  the  sky  wave  is  even  partially  horizontally  polarized. 

In  general,  sky  waves  do  have  horizontally  polarized  components. 

Thus  even  if  the  sky  wave  and  the  ground  wave  have  the  same  azimuthal 
angle  of  arrival  (which  is  not  always  true  either)  the  loop  will  not  give 
a  satisfactory  indication  this  angle. 

The  initial  solution  to  this  problem  lay  in  constructing  a  direction¬ 
finding  antenna  that  would  respond  only  to  the  vertically  polarized 

2 

component  of  the  wave.  The  Adcock  system,  as  it  is  called  after  its 
mventor,  does  just  this.  In  one  form  it  consists  of  two  vertical  elec¬ 
tric  dipoles  placed  parallel  to  each  other  so  that  they  and  the  support 
structure  form  the  letter  H.  The  dipoles  are  connected  180°  out  of 
phase  and  the  structure  is  rotated  about  the  center  of  the  cross  bar  of 
the  H.  If  the  spacing  between  the  dipoles  is  small  enough  (less  than 
1/4  wavelength)  the  receiving  pattern  has,  in  azimuth,  the  form  of  a 
rotating  figure  8.  The  D/F  is  performed  on  the  nulls  with  some 
auxiliary  antenna  to  remove  the  180°  ambiguity.  So  if  the  ground  wave 
and  sky  wave  have  the  same  azimuthal  angle  of  arrival  the  Adcock  is 
capable  of  determining  it. 

If  these  two  waves  do  not  have  the  same  azimuthal  angle  of  arrival, 
however,  the  bearing  indication  given  by  the  Adcock  may  have  a  large 
systematic  error.  In  fact  all  D/F  systems  which,  like  the  Adcock, 
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employ  directional  antennas  of  small  aperture,  are  subject  to  it.  It 

was  first  described  and  its  cause  explained  in  1923  by  Heiligtag  ^ 

4 

Recently  Hayden  reviewed  its  derivation.  A  summary  of  his  presen¬ 
tation  follows;  In  Figure  1-1  is  shown  the  standing  wave  field  produced 
by  two  vertically  polarized  plane  waves  whose  wave  normals  lie  in  the 
plane  of  the  paper  and  have  a  30°  angle  between  them.  The  equiphase 
surfaces  of  this  field  are  corrugated  (or,  if  the  two  waves  are  of  equal 
strength,  broken).  From  the  figure  it  can  be  seen  that  the  normal  to 
the  phase  front  varies  from  a  direction  between  the  wave  normals  (but 
closer  to  the  normal  to  the  stronger  wave),  to  a  direction  outside  the 
angle  between  the  wave  normals  (but  away  from  the  normal  to  the 
weaker  wave).  At  present  all  small  aperture  D/F  systems  give  an 
indicated  bearing  which  is  equal  to,  or  very  nearly  equal  to,  the  normal 
to  the  equiphase  surface.  Thus  the  indicated  bearing  of  such  a  system 
is  significantly  dependent  on  where  it  is  placed. 

One  answer  to  this  problem  is  to  construct  a  large  aperture  system 

which  will  average  out  the  "wrinkles"  in  the  wavefront  to  give  a  true 

bearing  indication.  In  principle,  this  is  what  the  Wullenweber^ 

6  7 

system,  the  cyclical -differential -phase  system  and  the  Doppler 
system  do.  Another  approach  is  to  average  over  the  bearings  of  a 
number  of  spatially  dispersed  small  aperture  systems.  Both  of  these 
methods  are  expensive  to  implement,  however,  so  considerable  effort 
has  gone  into  the  analysis  of  multipath  propagation,  with  the  hope  that 
cheaper  and  less  involved  D/F  systems  will  result. 

Because  multipath  propagation  in  the  ionosphere  is  the  major  source 

8-15 

of  error  in  the  HF  (1  -  30  MHz)  band,  considerable  effort  has  gone 
into  investigating  it  The  bulk  of  this  work  is  necessarily  experimental, 
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FIGURE  1-1  WAVE  INTERFERENCE  FIELD  FOR  RAY  SEPARATION  OF  3 CP 
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since  the  mechanisms  involved  are  quite  complicated.  Useful  results 
have  followed,  particularly  in  the  substitution  of  time  averaging  for 
space  averaging.  The  dynamic  nature  of  the  ionosphere  makes  this  a 
profitable  procedure.  Relatively  little  effort,  however,  has  gone  into 
the  second  source  of  multipath  propagation,  namely,  the  reflection  of 
electromagnetic  waves  from  surface  obstacles  in  the  vicinity  of  the 
direction-finder . 

Siting  error,  as  it  is  called,  has  been  overshadowed  by  the  larger 
errors  due  to  multipath  propagation  in  the  ionosphere,  but  it  is  by  no 
means  insignificant.  It  is  the  major  source  of  D/F  error  in  the  VHF 
and  UHF  bands  and  it  is  important  in  the  band  for  short  range 
(0  -  50  miles)  direction-finding.  The  work  to  follow  will  concentrate 
on  the  latter  area,  but  it  is  equally  applicable  to  the  VHF  and  UHF 
areas  as  well. 

1 6  25 

Previous  investigators  have  concentrated  on  the  effect  of  bearing 

accuracy  of  one  or  two  "point  reflectors"  (i.e.  ,  structures  whose 
major  dimension  is  much  smaller  than  a  wavelength,  or  thin  vertical 
towers)  in  the  vicinity  of  the  direction-finder.  The  comparatively 
simple  nature  of  the  reradiated  field  from  a  "point"  source  or  a  linear 
reflector  parallel  to  the  polarization  vector  leads  to  explicit  expressions 
for  the  expected  bearing  error  in  their  presence.  Needless  to  say, 
sites  where  one  or  two  dominant  point  reflectors  constitute  the  major 
source  of  reradiation  are  very  much  in  the  minority.  Far  more 
common  are  sites  characterized  by  several,  more  or  less,  extended 
sources  of  reradiation  (i.e.,  buildings,  hills,  clumps  of  trees).  Any 
object  whose  horizontal  dimensions  are  large  will  have  a  complex 
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induced  current  distribution  on  and  within  it.  The  reradiated  fields 
are  correspondingly  complex.  This  work  attempts  a  unified  analysis 
of  sites  possessing  obstacles  of  this  type  and,  by  adopting  a  statistical 
model,  arrives  at  a  single  parameter  useful  in  evaluating  them.  This 
parameter  can  be  measured  quickly  and  accurately  on  a  given  site, 
through  a  new  method  presented  here. 

The  statistical  techniques  given  have  not  been  applied  to  direction¬ 
finding  before.  Their  advantage  is  two-fold. 

1.  The  problem  of  site  evaluation  is  given  a  more  logical 
solution,  in  place  of  the  rule-of-thumb  methods  in  use 
at  present. 

2.  The  preliminary  measurement  procedure  necessary  for 
site  selection  is  considerably  simplified. 

It  is  hoped  that  the  conclusions  presented  will  provide  a  practical 
guide  for  engineers  engaged  in  choosing  D/F  sites. 
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Chapter  2 

DEVELOPMENT  OF  AN  ANALYTICAL  MODEL  OF  A  D/F  SITE 

A.  Hypothesis  to  Explain  D/F  Errors 

About  twenty  years  ago,  various  investigators,  encouraged  by  their 
success  with  shipboard  direction-finders,  were  examining  the  pos¬ 
sibility  of  "calibrating"  a  D/F  site  by  plotting  measured  errors  in 
bearing  against  ang^e  of  arrival  (for  a  fixed  frequency).  One  of  these 
investigators,  W.  Ross,  *  performed  a  series  of  experiments  on 
several  sites  to  determine  if  calibration  were  feasible.  The  sites  he 
selected  were,  in  general,  very  good,  as  far  as  could  be  judged  by 
visual  inspection;  i.e.,  the  direction-finding  stations  were  situated 
in  the  middle  of  fairly  flat  terrain.  Small  hills  and  folds  in  the 
ground  were  present,  but,  from  past  experience,  did  not  appear  to 
be  sufficient  reason  to  reject  the  site.  The  procedure  he  followed 
was  to  examine  the  size  of  the  bearing  error  as  a  function  of  each  of 
three  variables,  namely,  azimuthal  position  of  a  target  transmitter, 
frequency,  and  distance  between  the  direction-finder  and  the  target 
transmitter.  The  experiments  were  carried  out  at  various  frequen¬ 
cies  between  6  and  15  MHz. 


Ross  found  that  site  calibration  was  not  feasible.  While  it  was  poss¬ 
ible  to  measure  bearing  error  as  each  of  the  above-mentioned 
quantities  was  varied  in  turn,  the  amount  of  information  required 
for  a  complete  calibration  was  too  great  to  obtain  in  a  reasonable 
period  of  time.  This  conclusion  followed  from  three  measured 
effects : 
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1.  The  observed  errors  were  not  predictable  functions  of  the 
azimuthal  position  of  the  target-transmitter. 

2.  Calibration  at  one  frequency  gave  no  information  about 
calibration  at  another.  There  were,  for  example,  signifi¬ 
cant  differences  in  observed  errors  when  the  frequency  was 
changed  by  only  two  per  cent. 

3.  There  was  a  lesser,  but  noticeable  dependence  of  bearing 
error  on  the  distance  between  the  target  transmitter  and  the 
direction-finder.  This  dependence  is  less  strong  as  the 
distance  increases. 

The  erratic  character  of  the  bearing  errors  measured  led  Ross  to 
hypothesize  that  they  were  random  in  character;  i.  e.  ,  "for  any  par¬ 
ticular  bearing  and  wavelength  the  error  which  might  be  observed 
could  be  regarded  as  the  total  effect  of  a  large  number  of  causes  each 
contributing  a  small  fraction  to  the  total  error,  the  causes  of  the 
errors  being  regarded  as  scattered  at  random  over  a  considerable 
area  all  around  the  direction-finder."  He  further  suggested  that  "in 
addition  to  the  direct  wave  from  the  transmitter,  a  large  number  of 
reflected  waves  from  many  directions  originate  from  countless  suit¬ 
able  reflectors;  e.g.,  trees,  metal  poles  and  pylons,  buildings  with 
metal  frameworks,  small  hills,  etc.,  scattered  over  a  wide  area 
surrounding  the  direction-finder." 

Ross  presented  some  experimental  justification  for  his  hypothesis; 
e.g.  ,  it  is  a  good  explanation  of  the  very  variable  nature  of  the  error 
with  small  changes  in  wavelength.  He  did  not,  however,  attempt  to 
express  his  ideas  mathematically  and  to  derive  therefrom  the  observ- 
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able  results.  The  modern  statistical  approach  which  follows  is  an 
application  of  the  theory  of  rough-surface  scattering  to  the  D/F 
problem. 


In  subsequent  sections  of  this  chapter  a  statistical  model  of  a  site  is 

2  3 

developed  based  on  the  work  of  Beckmann  and  Siddiqui  and  a  single 
parameter  is  suggested  as  a  means  of  evaluating  a  site.  The  proba¬ 
bility  of  errors  in  bearing  of  a  specified  magnitude  depends  both  on 
this  parameter  and  on  the  type  of  direction-finder  used.  In  Appen¬ 
dix  A  the  relation  of  this  parameter  to  bearing  uncertainty  for  an 
Adcock  type  direction-finder  is  presented  A  similar  relationship  may 
be  worked  out  for  other  types  of  D/F  systems. 

B.  General  Assumptions 

Certain  general  assumptions  are  necessary  in  order  to  develop  the 
analysis.  Before  stating  these,  however,  the  problem  is  defined 
more  precisely  by  the  following  observations: 

1.  Ross  examined  the  influence  of  three  variables  on  bearing 
error.  The  first  area  considered  in  the  analysis  is  the  in¬ 
fluence  of  only  one  of  these  variables,  angle,  on  bearing 
error.  The  necessary  assumptions  in  this  case  seem  more 
reasonable  (i.e.  ,  defensible)  than  in  the  case  of  the  other 
two  variables . 

2.  The  emphasis  is  placed  upon  the  relationships  between  the 
reradiated  and  the  direct  waves,  without  reference  to  any 
particular  type  of  direction-finder.  This  is  done  in  order 
to  make  the  results  as  general  as  possible.  The  effect  of 
the  site  on  the  spatial  structure  of  the  fields  is  independent 
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of  the  direction-finder,  but  the  bearing  errors  resulting 
from  site  imperfections  are  different  for  different 
direction-finders . 

3.  The  primary  concern  is  with  sites  which  are  characterized 
by  the  presence  of  extended  reradiators.  Most  sites  are  of 
this  type. 

The  assumptions  essential  to  the  analysis  may  now  be  stated.  These 
assumptions,  in  mathematical  terms,  are  those  employed  in  the 
theory  of  electromagnetic  scattering  from  a  rough  surface,  and  those 
used  in  the  theory  of  ionospheric  propagation  Their  interpretation 
in  view  of  the  present  problem  requires  some  words  of  justification. 
Each  listed  assumption  is  therefore  followed  by  an  explanation  sug¬ 
gesting  the  reason  why  it  may  be  applied. 

1.  The  site  is  illuminated  by  a  vertically  polarized  monochro¬ 
matic  plane  wave. 

This  assumption  restricts  the  analysis  to  ground  waves 
launched  by  a  vertically  polarized  transmitting  antenna 
placed  on  the  surface  of  the  earth.  Strictly  speaking,  of 
course,  there  is  always  present  in  a  ground  wave  a  small 
horizontally  polarized  component,  whose  magnitude  depends 
on  the  conductivity  of  the  earth.  A  vertically  polarized  re¬ 
ceiving  antenna,  however,  will  not  respond  to  this  component 
and,  in  effect,  will  behave  as  if  the  wave  were  vertically 
polarized.  This  is  not  to  say  that  the  methods  to  be  employed 
are  not  applicable  to  other  polarizations.  In  the  conclusions 
it  is  suggested  how  this  might  be  accomplished. 


2  For  any  azimuthal  position  of  the  transmitting  antenna,  the 
field  at  the  direction-finder  is  the  resultant  field  of  a  direct 
wave  from  the  transmitter  and  a  number  of  scattered  waves. 

The  scattered  waves  refer  to  the  fields  radiated  by  the  sur¬ 
face  and  volume  currents  induced  by  the  direct  wave  on 
various  obstacles  present  on  the  D/F  site, 

3.  At  the  direction-finder  the  amplitude  and  phase  (referred  to 
the  direct  wave)  of  each  of  the  scattered  waves  varies  in  a 
random  manner  as  the  angle  of  arrival  of  the  direct  wave 
changes 

The  induced  current  distribution  on  an  irregularly  shaped 
obstacle,  whose  horizontal  dimensions  are  large,  has  a 
complicated  form.  As  the  illumination  angle  of  the  direct 
wave  changes,  this  current  distribution  changes  also. 

Given  an  obstacle,  it  is  a  formidable  job  to  calculate  this 

distribution  for  a  specific  angle  of  arrival.  It  is  virtually 

,  4 

impossible  to  do  so  for  all  angles.  Cramer  points  out  that 

in  such  a  situation  a  probabilistic  model  may  be  more  useful 

than  a  deterministic  model.  He  states  that  the  observable 

characteristics  of  a  physical  phenomenon  may  exhibit  the 

fluctuations  typical  of  random  variables,  even  when  the  laws 

of  the  phenomenon  are  regarded  as  well  known,  provided 

that  these  laws  are  sufficiently  complicated.  He  gives  the 

5 

following  example: 

"We  do  assume  that  it  is  possible  to  predict  the  annual  num¬ 
ber  of  eclipses  (of  the  sun),  and  if  the  requisite  tables  are 
available,  anybody  can  undertake  to  make  such  predictions. 
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Without  the  tables,  however,  it  would  be  rather  a  formidable 
task  to  work  out  the  necessary  calculations,  and  if  these  dif¬ 
ficulties  should  be  considered  insurmountable,  prediction 
would  still  be  practically  impossible,  and  the  fluctuations  in 
the  annual  number  of  eclipses  would  seem  comparable  to  the 
fluctuations  in  a  sequence  of  games  of  chance.  " 

It  is  convenient,  then,  to  consider  the  amplitude  and  phase 
of  the  scattered  wave,  which  are  observable  quantities  pro¬ 
duced  by  the  induced  current  distribution  on  an  obstacle,  as 
random  variables.  It  is  difficult  to  say  when  such  a  model 
may  not  apply.  It  does  not  seem  reasonable,  for  example, 
when  the  obstacle  is  a  thin  tower  of  known  height.  Then  the 
amplitude  and  phase  of  the  scattered  waves  appear  to  be 
more  or  less  predictable  functions  of  the  angle  of  arrival. 

The  question  of  the  applicability  of  the  model  in  this  and 
other  cases  is  dealt  with  more  fully  in  subsequent  sections 
of  this  work  and  particularly  in  the  CONCLUSIONS  (Chapter  5). 

4.  The  number  of  scattered  waves  is  sufficiently  large. 

How  many  waves  will  be  considered  sufficient  depends  to  a 
large  extent  on  the  probability  density  distributions  of  the 
amplitudes  and  phases  of  the  scattered  waves.  In  certain 
cases  only  one  may  be  required.  In  most  practical  cases, 
no  more  than  six^  are  required.  Further  explanation  will 
be  given  in  the  next  section. 

These  four  assumptions  are  sufficient  to  derive  the  probability 
density  distribution  of  the  amplitude  of  the  field  at  the  direction¬ 
finder.  This  distribution,  however,  is  a  complicated  expression, 
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depending  on  four  parameters  which  are  not  useful  quantities  for 
site  evaluation  In  order  to  derive  practical  results  from  the  analy¬ 
sis  ,  three  more  assumptions  will  be  made  The  first  two  have  good 
justification  The  last  is  less  justifiable,  but  an  indication  of  how  it 
may  be  generalized  and  therefore  better  justified,  will  be  given. 

5  Each  of  the  scattered  waves  is  statistically  independent  of 
the  others 

This  assumption  implies  that  the  obstacle -to-obstacle  inter¬ 
action  effects  are  small  compared  to  the  obstacle-to-direct 
wave  interaction  effect  In  other  words,  the  current  induced 
on  an  obstacle  A  by  the  current  induced  on  obstacle  B  is 
small  compared  to  the  current  induced  on  obstacle  A  by  the 
direct  wave 

6  The  amplitude  and  phase  of  each  scattered  wave  are  statis¬ 
tically  independent  random  variables. 

From  assumption  3  it  can  be  argued  that  amplitude  and 
phase  may  be  considered  to  be  independent  A  change  in  the 
induced  current  distribution  on  an  obstacle  will  not  affect 
the  relationship  between  them  in  an  easily  predictable  way 

7  The  phase  of  any  one  of  the  scattered  waves  is  uniformly 
distributed  between  -tt  and  tt. 

This  assumption  gives  rise  to  a  particularly  useful  result. 

A  more  general  assumption  would  be  that  the  probability 
density  distribution  of  the  phase  is  symmetrically  distri¬ 
buted  about  zero  This  case  will  be  fully  examined  in  the 
next  section  It  will  be  shown  however,  that  the  difference 


between  the  result  of  assumption  7,  and  the  result  of  assuming 
a  symmetric  distribution  is  not,  in  general,  significant  for 
the  actual  range  of  parameters  encountered. 


C.  Probability  Distribution  of  Electric  Field  Amplitude 

The  vertically  polarized  field  at  the  direction-finder  is  the  sum  of  the 

vertically  polarized  components  of  the  individual  waves  scattered  by 

various  obstacles  and  the  direct  wave.  If  the  scattered  waves  are 

j$ 

represented  by  complex  numbers  E  e  n  (suppressing  the  time 

t  n 

factor  e^  ,  as  usual),  where  the  values  of  $  are  referred  to  the 

n 

direct  wave,  each  individual  wave  will  be  represented  by  a  vector  in 
the  complex  plane.  The  resultant  field  will  be  given  a  vector  sum  in 
the  complex  plane  (Figure  2-1)  This  sum  is  expressed  mathemati¬ 
cally 

N 
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(Eq  2-1) 


where  E  is  the  amplitude  of  the  direct  wave  and  it  is  understood  that 
d 

the  quantities  E^,  $  ,  the  amplitude  and  phase  of  the  nth  scattered 

wave  are  random  variables  whose  magnitude  fluctuates  as  the  azi¬ 
muthal  angle  of  arrival  of  the  direct  wave  changes  Under  the  as¬ 
sumptions  given  in  the  preceding  section,  the  probability  density 
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distribution  of  E  is  a  Rice  distribution,  given  by 
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where  <E  >  is  the  sum  of  the  mean  square  amplitudes  of  the  scattered 
s 

waves,  averaged  over  all  possible  values  of  angle,  i.  e.  : 


N  N  N 

2  r-»  O  n  V  2 

<E  >  =  <(  V  £  e  )  {  )  E  e  )>  =  )  <E  >  (Eq.  2-3) 
s  L,  n  Z_i  n  L,  n 
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(since  p( I  )  =  p( $  )  =  1/2tt  for  all  n,  m)  and  where  I  is  the  modified 

n  m  o 

Bessel  Function  of  the  first  kind  and  Oth  order. 


As  was  pointed  out  in  the  previous  section,  the  comparatively  simple 

form  of  the  distribution  is  due  primarily  to  assumption  7.  Adopting 

the  more  general  assumption  that  the  p( i  )  are  all  symmetrical  about 

8  n 

i  =0,  Beckmann  has  derived  an  expression  for  p(E).  A  summary 
n 

of  his  derivation  follows. 


First,  four  quantities  must  be  defined: 

1)  The  average  value  of  the  real  part  of  E  e 


i* 


N 

a  =  <E  cos  $>  =  E  +  <  ^  E  cos  $  > 
d  L  n  n 
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(Eq.  2-4) 


2)  The  average  value  of  the  imaginary  part  of  E  e 


j* 


N 


3  =  <E  sin  $>  =  <  '  E  sin  $  > 
L,  n  n 


(Eq  .2-5) 


n  =  1 


3)  The  variance  of  the  real  part  of  E  e' 


i# 
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N  N 

2  2  ^  r-,  2 

s ,  =  < (E  cos  $)  >  -  a  =  <  (  E  cos  $)>-<)  E  cos  $  > 

1  L,  n  n  u  n  n 

n=l  n=l 

(Eq.  2-6) 

4)  the  variance  of  the  imaginary  part  of  E  eJ 

N  N 

2  2c  2  c  2 

s_  =  <(E  sin  $)  >  -  0  =  <(  )  E  sin  $)>-<)  E  sin  $  > 

2  L>  n  n  L  n  n 

n-1  n=l 

(Eq.  2-7) 


Now  if  N,  the  number  of  scattered  waves  is  large  enough  and  the  dis¬ 
tributions  of  E  cos  $  and  E  sin  $  satisfy  the  conditions  of  the 
n  n  n  n 

Central  Limit  Theorem,  then  E  cos  $  and  E  sin  $  will  be  normally 
distributed  with  mean  values  a  and  0  and  variances  s^  and  s^  re¬ 
spectively.  Moreover,  if  the  distributions  p(  $  )  are  symmetrical 

n 

about  zero  as  we  have  assumed,  then  the  two  quantities  E  cos  $  and 
E  sin  $>  will  be  uncorrelated.  The  proof  of  this  proceeds  as  follows: 

The  covariance  of  these  quantities  is  defined  as 

C  =  <(E  cos  $)  (E  sin  $)>  -  0.0 

N  N 

=  <(E  +  V  E  cos  $  )  (  ^  E  sin  $  )>  -  00  (Eq.  2-8) 

dZ_in  nZ_,n  n 

n= l  n= 1 


But  for  p($n)  symmetrical  about  zero 
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1^1 

G  =  <)  E  sin  $  >  =  0 
Li  n  n 


hence 


C  =  < 


N  N 

I  5>. 

n=l  m=l 


E  cos  i  sin  $  > 
m  n  m 


If  the  identity  is  employed  that 


cos  i  sin  $ 
n  m 


(sin  ($  +  $  )  -  sin  ($  -  $  )  ) 

2  n  m  n  m 


it  follows  that 


C  =  0 


where  the  assumption  that  $  , 
been  used. 


(Eq.  2-9) 


.  .  ,  are  independent  variables  has 


Since  E  cos  $  and  E  sin  $  are  each  normally  distributed  and  uncor¬ 
related,  they  are  statistically  independent.  Thus  the  two-dimensional 
distribution  of  E  cos  §  and  E  sin  i  is  the  product  of  their  individual 
distributions : 


p(E  cos  $  E  sin  $)  =  -  —  exp 

2tt  7  s  js2 


(E  cos  $  -  a 
2s . 


(E  sin  $) 
2s  „ 


(Eq.  2-10) 


In  order  to  determine  p(E),  the  distribution  function  of  the  amplitude 
of  the  resultant,  the  variables  in  (2-10)  must  be  changed  from 


(E  cos  $,  E  sin  $)  to  (E,  $).  The  Jacobian  of  this  transformation  is  E 
hence  the  distribution  p(E)  is  given  by 


tt  n 

p(E)  =  ^  p(E,  §)  d§  =  E  C  p(E  cos  $,  E  sin  $)  d$ 
-n  -rr 


TT 


-TT 


(E  cos  $  -  q) 
2S! 


(E  sin  j) 

2s„ 


d$ 


(Eq.  2-11) 


This  integral  is  solved  by  Beckmann  .  The  final  form  is 


(Eq.  2-12) 


where  I  is  the  modified  Bessel  Function  of  order  m  and 
m 


e  =  1  rr.  =  0 

rr 

=  2  m  i  0  (Eq.  2-13) 


If  assumption  7  is  introduced,  (2-12)  becomes  a  Rice  distribution,  for 
then 
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(Eq.  2-14) 


giving 


N 


_2  1  _2 

< E  >  = - <E  > 
n  2  s 


n=  1 


N 
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2  1  2 

<E  >  =4< E  > 
n  2  s 


1  t-2 

=  T  <  E  > 
2  s 


n=  1 


p(E)  = 


2E 


<EZ> 

s 


exp 


2  2 
'<Ed  +  E 


<e2> 

s 


4  /  2Ea 

4  (<e?> 


(Eq.  2-15) 


the  distribution  stated  at  the  beginning  of  this  section. 


Section  B  mentioned  that  for  the  range  of  parameters  under  consid¬ 
eration  in  the  present  problem,  there  is  little  difference  between 
assuming  that  the  p(  $  )  are  symmetrical  about  zero  and  assuming 

that  the  p($  )  =  1/2tt  from  -tt  to  tt  .  It  is  expected  that  the  p(§  )  will 
n  n 

not  be  unusual  (e.g.  ,  5  -function)  distributions,  so  s^  and  s^  are 

probably  not  greatly  different  from  each  other.  In  this  case  the 

dominant  term  in  the  summation  of  (2-12)  will  be  the  first,  I  (aE/s,) 

o  1 

since 


/  1 

!m  (  4s  j 


-  s 


-V 


is  small  for  m>0.  The  result  is  that  (2-12)  and  (2-15)  have  approx¬ 
imately  the  same  form  regardless  of  which  of  the  two  assumptions 
are  made. 
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An  observation  should  be  made  on  how  large  N  must  be  to  satisfy  the 

conditions  of  the  Central  Limit  Theorem.  This  depends  on  the  density 

distributions  of  the  quantities  being  summed.  If  E  cos  $  and  E  sin 

n  n  n 

$  are  distributed  normally,  E  cos  $  and  E  sin  $  are  normally  dis¬ 
tributed,  then,  for  N  2  1.  If  they  are  distributed  uniformly,  E  cos  $ 
and  E  sin  $  are  for  all  practical  purposes  distributed  normally  for 

N  a  4.  Even  if  the  E  were  constant,  with  $  a  random  variable  uni- 

n  n 

formly  distributed  between  -it  and  tt  ,  then  it  has  been  shown  that 
E  cos  $  and  E  sin  $  are  normally  distributed  for  N  £  6.  ^ 


An  approximation  of  (2-15)  useful  for  further  work  is  now  made.  First, 
the  expression  is  modified  by  carrying  out  the  substitutions 


r  = 


(Eq.  2-16) 


This  gives 


p(r) 


2r 

—  exp 


U  *  rZ) 

k2 


(Eq.  2-17) 


It  is  expected  that  the  magnitude  of  the  scattered  waves  will  be  small 

compared  to  the  magnitude  of  the  direct  wave.  Mathematically  this 
2 

means  that  k  <  <  1  and  r  «  1 ,  so  it  is  useful  to  find  an  approx¬ 
imation  for  (2-17)  in  this  range.  For  large  values  of  the  argument 
the  Bessel  Function 


(Eq.  2-18) 
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Substituting  (2-18)  into  (2-17), 


p(r)  *>  ■  —  exp 

A? 


(r  -  I)' 


(Eq.  2-19) 


is  obtained.  Therefore  for  small  values  of  k  and  r  «  1,  the  dis¬ 
tribution  function  is  approximately  normally  distributed  with  a  mean 

2  . 

of  1  and  a  variance  of  k  / 2 . 

2 

It  is  suggested  that  k  be  used  as  the  evaluation  parameter  for  sites 

2 

which  meet  the  assumptions.  Values  of  k  above  an  arbitrary  cri¬ 
terion  will  indicate  an  intolerable  level  of  scattered  power  present  on 
the  site.  The  utility  of  this  parameter  may  be  gauged  from  the  dis¬ 
cussion  of  narrow-aperture  direction-finders  given  in  Appendix  A. 


D.  Summary 

In  this  chapter  the  theory  of  rough-surface  scattering  is  employed  to 
develop  an  analytical  model  which  will  be  as  general  as  possible  in  its 
application  to  D/F  sites.  A  hypothesis  based  on  measurements  per¬ 
formed  on  several  sites  is  presented  in  Section  A.  In  Section  B,  part 
of  this  hypothesis  is  quantized  by  proposing  in  mathematical  terms,  a 
number  of  general  assumptions  concerning  the  nature  of  the  reradiated 
fields.  In  Section  C  the  consequences  of  these  assumptions  are  de¬ 
rived  and  a  single  parameter,  which  is  the  ratio  of  scattered  power  to 
direct  power,  to  characterize  D/F  sites  is  suggested. 
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Chapter  3 


SITE  EVALUATION  BY  FREQUENCY  VARIATION 

A  Scatter i n g  Characteristics  as  a  Function  of  Frequency 

In  the  statistical  model  of  rough- surfac e  scattering  (or  ionospheric 
propagation)  which  has  been  applied  to  the  problem  of  D/F  site  evalu¬ 
ation,  an  interpretation  is  given  of  the  variation  in  the  intensity  of  the 
total  electric  field  at  the  direction-finder  as  the  angle -of-arrival  of 
the  direct  wave  changes.  The  statistical  approach  to  this  problem, 
however,  leads  to  an  inference  not  discussed  with  relation  to  the 
rough-surface  scattering  problem  namely  that  the  fluctuations  in 
the  observed  electric  field  intensity  when  the  frequency  of  the  direct 
wave  is  varied  might  possess  the  same  statistical  characteristics  as 
when  the  angle-of-arrival  is  varied  This  implication  has  not  been 
drawn  in  the  rough- surface  scattering  model  or  the  ionospheric  propa 
gation  model,  since  both  are  primarily  concerned  with  time  or  angle 
varying  signal  fluctuations  at  a  fixed  frequency 

Clearly,  two  questions  are  being  posed: 

1.  Is  the  electric  field  intensity  at  the  direction-finder  a  ran¬ 
dom  function  of  both  angle-of-arrival  and  frequency? 

2.  Does  the  probability  distribution  which  describes  the 
variation- with-angle  case  also  describe  the  variation-with- 
frequency  case? 

Should  the  answer  to  both  questions  be  in  the  affirmative,  then  an  im¬ 
portant  result  follows;  i.e.  ,  that  the  suggested  site  evaluation  para- 
2 

meter,  k  ,  can  be  determined  by  examining  the  fluctuations  of  the 
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received  intensity  as  the  frequency  of  the  direct  wave  is  swept  over 
some  arbitrary  bandwidth,  the  angle -of-arrival  being  fixed. 

The  questions  above  are  addressed  in  this  chapter  in  two  ways. 

First  of  all,  in  this  section,  a  discussion  is  given  of  the  reasonable¬ 
ness  of  the  assumptions  of  Chapter  2,  as  they  apply  to  the  scattering 
behaviour  of  complex  obstacles  as  a  function  of  frequency.  Secondly, 
in  the  next  section,  the  answers  to  these  questions  are  given  for  sev¬ 
eral  arrays  of  hemispherical  bosses  on  a  conducting  ground  plane. 

Assumptions  1,  2,  and  4  of  Chapter  2  require  no  explanation  for  this 
case  of  variable  frequency;  the  others  do. 

1.  The  site  is  illuminated  by  a  vertically  polarized  monochro¬ 
matic  plane  wave. 

2.  For  any  transmitter  frequency,  the  field  at  the  direction¬ 
finder  is  the  resultant  field  of  a  direct  wave  and  a  number  of 
scattered  waves. 

3.  At  the  direction-finder  the  amplitude  and  phase  (referred  to 
the  direct  wave)  of  each  of  the  scattered  waves  varies  in  a 
random  manner  as  the  frequency  of  the  direct  wave  changes. 

The  same  justification  is  given  as  in  Chapter  2.  If  an  obstacle  is  suf¬ 
ficiently  large  with  respect  to  the  longest  wavelength  used,  then  the 
current  distribution  on  the  obstacle  changes  in  a  complicated  fashion 
as  the  frequency  of  the  direct  wave  changes.  According  to  Cramer, 
it  is  quite  reasonable  to  consider  the  amplitude  and  phase  of  the  field 
produced  by  this  complicated  current  distribution  as  random  variables 
dependent  on  frequency. 


32 


4.  The  number  of  scattered  waves  is  sufficiently  large. 

5.  Each  of  the  scattered  waves  is  statistically  independent  of 
the  others. 

If  the  obstacle -to-obstacle  distance  is  large  enough  with  respect  to 
the  longest  wavelength  used,  the  effects  of  multiple  scattering  may  be 
neglected,  or,  in  other  words,  the  induced  curre  it  on  one  obstacle 
has  negligible  effect  on  the  induced  current  on  another  obstacle.  It  is 
reasonable  to  assume  the  scattered  waves  are  independent  under  this 
condition. 

6.  The  amplitude  and  phase  of  each  scattered  wave  are  statisti¬ 
cally  independent  random  variables  over  frequency. 

For  complex  obstacles  the  amplitude  and  phase  of  each  scat¬ 
tered  wave  are  not  predictably  related  when  frequency  is 
varied.  From  Cramer's  point  of  view,  therefore,  if  it  is 
convenient  to  assume  them  independent  of  each  other,  one 
may  do  so;  the  assumption  does  not  contradict  any  observable 
behaviour. 

7.  The  phase  of  any  one  of  the  scattered  waves  is  uniformly 
distributed  between  -tt  and  IT. 

This  is  again  the  least  defensible  assumption.  As  before,  it 
can  be  generalized.  In  any  case,  it  is  reasonable  to  expect 
the  phase  of  a  wave  reradiated  from  an  obstacle  to  fluctuate 
over  a  period  as  the  frequency  changes,  as  long  as  the  range 
of  frequencies  is  large  enough. 
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Now  at  any  frequency,  the  scattered  waves  may  be  represented  by 

j$ 

complex  numbers  E  eJ  n  where  i  is  the  phase  referred  to  the 

n  n 

direct  wave.  The  sum  of  the  direct  wave  plus  the  scattered  waves  is 
expressed,  as  before, 


Eej$  =  Ed  +  £  E^e  n  (Eq.  3-1) 

n=  1 

where  E  is  the  amplitude  of  the  direct  wave  and  E  i  ,  the  am- 
d  n  n 

plitude  and  phase  of  the  nth  scattered  wave  are  random  variables 
whose  magnitude  fluctuates  as  the  frequency  of  the  direct  wave 
changes.  Under  the  assumptions  given  above,  the  random  fluctua¬ 
tions  of  E  as  frequency  changes  are  described  by  the  Rice  distribu¬ 
tion: 


P(E) 


2E 

- ; -  exp 

<E  > 
s  0) 


(E2  +  E2) 
<e2> 

s  uu 


I 

o 


2E 


<e2> 

s  (JU 


1 


(Eq.  3-2) 


where  <E  > 


s  uu 

tered  waves, 


is  the  sum  of  the  mean  square  amplitudes  of  the  scat- 
averaged  over  frequency;  i.  e.  , 


N 

<e2>  =  y  <e2> 

s  uu  l—*  n  uu 

n=  1 


(Eq.  3-3) 
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If  the  substitutions 


#-  ;  K2 

Ed 


_  2 

<E  >  uj 
s 


(Eq.  3-4) 


are  made,  then 


/  \  2r 
p(r)  =  —r  exp 

k 

OJ 


|1  4  r 
.  2 


1  H 

°  k2 

UJ 


(Eq.  3-5) 


The  assumptions  lead  to  the  same  type  of  distribution  as  that  given  in 

2 

Chapter  2,  with  the  difference  that  the  decribing  parameter  k  of 

2  W 

(3-5)  is  not  necessarily  the  same  as  the  k  of  (2-16).  Referring  to 

2  2  2  2 
the  k  of  (2-16)  as  k  (angle)  and  the  k^  of  (3-5)  as  k  (frequency), 

the  inference  made  at  the  beginning  of  this  section  is  that 


2  2 
k  (angle)  =  k  (frequency) 


(Eq.  3-6) 


Or,  in  other  words,  the  two-dimensional  variation  of  r  (with  angle  and 
frequency)  can  be  described  by  the  same  distribution. 

For  complex  obstacles  it  is  a  formidable  job  to  investigate  analyti¬ 
cally  the  conditions  under  which  the  inference  might  be  true.  If  it 
can  be  shown,  however,  that  the  inference  is  applicable  to  obstacles 
of  simple  shape,  then  it  may  be  applied  with  more  confidence  to  com¬ 
plex  obstacles.  In  the  next  section  its  applicability  to  a  randomly 
distributed  array  of  perfectly  conducting  hemispherical  bosses  on  a 
perfectly  conducting  ground  plane  is  examined  in  detail. 
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B.  Scattering  from  Hemispherical  Bosses  on  a  Ground  Plane 

As  was  shown  by  Beckmann^  scattering  from  a  perfectly  conducting 
hemispherical  boss  on  a  ground  plane  is  equivalent  to  the  scattering 
from  a  sphere  of  the  direct  wave  and  its  image  in  the  plane.  Accord¬ 
ingly,  the  specific  case  to  be  examined  is  a  planar  array  of  randomly 
distributed  perfectly  conducting  spheres  of  different  radii.  A  plane 
wave,  whose  wave  normal  is  parallel  to  the  plane  of  the  array,  and 
whose  electric  field  vector  is  perpendicular  to  the  plane  of  the  array, 
impinges  upon  it.  The  situation  is  illustrated  in  Figure  3-1. 


If  the  spheres  are  placed  sufficiently  far  apart  so  that  the  effects  of 
multiple  scattering  can  be  neglected,  the  total  scattered  electric 
field  at  some  point  (P  in  3-7)  within  the  array  is  the  sum  of  the  scat¬ 
tered  fields  from  each  of  the  spheres  considered  individually,  with 
a  suitable  adjustment  in  phase.  Referring  to  the  figure,  the  scattered 
field  at  point  P  of  the  ith  sphere  of  radius  a.  located  at  p.,  i|i.,  with 

respect  to  the  x  -  y  coordinate  system  is  denoted 

E.  (0.,  p.,a.) 

where  0^  is  the  smaller  angle  between  the  wave  normal  and  the  radius 
vector  from  P  to  the  center  of  the  sphere.  If  iji  is  the  angle  between 
the  wave  normal  and  the  x-axis,  then 


0.  =  | TT  +  (ilf  -  i|/.)  |  if  ilf  <  iji. 

=  |tt  -  (t  -  *.)  |  if  ♦  > 


(Eq.  3-7) 


(0.  is  restricted  to  the  range  0  -  TT  for  computational  purposes.) 
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WAVE 


FIGURE  3-1 


A  Z-POLARIZED  PLANE  WAVE  OF  UNIT  AMPLITUJ 
INCIDENT  ON  A  PLANAR  ARRAY  OF  PERFECTLY  ( 
DUCTING  SPHERES.  (E  POINTS  OUT  OF  THE  PAGE 


The  total  scattered  field  at  P  is  thus  given  by 
N 

E  =  Y  exp  [ - jk o  .  cos  (i<  -  ’tj.)  ]  E.  (9.  ,p.,a.)  (Eq.  3-8) 
t  Z_,  i  11111 

i=  1 

where  N  is  the  total  number  of  spheres,  and  ft  -  2n/\  is  the  wave 
number  of  free  space. 


~  2,3 

An  expression  for  E.(0^,  p.,a.)  has  been  given  by  several  investi¬ 
gators.  If  the  magnitude  of  the  incident  field  is  unity,  then 


i  '  dm  d  pl(c°s  0.) 

r  1  V  i  u'2'  v  •  o  n  i 

E.  =  -  - -  >  I  c  H  (ft  o.)  sin  9. 

i  np.Z_.inn 


1  n  =  1  L 


i  d(cos  0.) 


..  d  A(2) 

-  jb  T7T~T  H  (ftb  ) 
n  d( n  o. )  n  i 


P  (cos  9.) 
n _ i 

sin  0. 

l 


(Eq.  3-9) 


where 


.  -n. 


J  (ft  a.) 

n  i 


n(n  +  1) 

AU)  i 

H  (ft  a.) 

n  l 

j  n(2n  +  1) 

d  A 

.  J  (ft  a.) 

d(u.)  n  i 

i 

n(n  +  1) 

d  e,<z> 

d(fta.)  n  '  T 

l 

A  A(2) 

J  (x),H  (x)  are  spherical  Bessel  Functions  of  order  n 


P^(cos  0)  are  Legendre  Functions  of  the  first  kind  of  order  n. 
n 
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A  computer  program  was  written  to  calculate  E^.  As  a  check  on  the 

program,  the  normalized  bistatic  scattering  cross  sections  of  a  single 

sphere  were  computed  for  40  values  of  ft  a.  equispaced  from  ft  a.  =  .  65 

tofta.  =  6.5.  The  bistatic  scattering  cross  section  a  is  defined  as 
i  P 

a  (8.)  =  lim  4n  p2  I  E.  (9.)  I  2 
p  i  pr*00  i  ii 

It  is  the  ratio  of  the  total  power  reradiated  by  a  fictitious  isotropic 
scatterer  (that  maintains  the  same  field  in  all  directions  as  that  main¬ 
tained  by  the  sphere  in  a  specified  direction)  to  the  real  magnitude  of 
the  Poynting  vector  of  the  incident  wave  at  the  sphere. 

A  sample  result  of  the  calculations  is  given  in  Figure  3-2.  The 

values  of  fta^  =  1.  7,  and  2.  3  are  those  for  which  the  value  of  the  back 

scattering  cross  section  (o(18Cf  )  )  is  a  maximum  or  a  minimum  when 

plotted  as  a  function  of  fta^.  The  curves  calculated  for  all  maximum 

and  minimum  values  of  a (180°)  are  the  same  as  those  given  by  King 
4 

and  Wu.  Since  the  bistatic  scattering  cross  section  is  symmetric 

in  8  it  is  presented  only  for  0 .  =  0  to  180° . 
i  i 

The  values  of  ft  a^  for  which  the  bistatic  scattering  cross  sections 
were  calculated  lie  in  the  "resonance"  region  of  the  sphere,  charac¬ 
terized  by  significant  fluctuations  in  the  amplitude  and  phase  of  the 
scattered  field  as  the  angle  of  arrival  and  the  frequency  of  the  direct 
wave  changes.  Although  the  behaviour  of  the  scattered  wave  in  this 
region  may  not  completely  satisfy  the  assumptions  made  in  Chapters 
2-B  and  3-A,  the  resonance  region  is  the  one  likely  to  give  the  best 
results  in  a  statistical  analysis.  Accordingly,  all  computations  were 
carried  out  in  this  range  of  fta^. 
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FIGURE  3-2  NORMALIZED  BISTATIC  SCATTERING  CROSS  SECTIONS  OF  A 

PERFECTLY  CONDUCTING  SPHERE  IN  A  PLANE-WAVE  FIELD  *■ 
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In  Figures  3-4,  3-6,  3-8,  3-10,  3-12  and  3-14,  are  given  plots  of  the 
total  (direct  plus  scattered)  electric  field  amplitude  as  the  angle  of 
arrival  (constant  frequency)  and  the  frequency  (constant  angle  of  ar¬ 
rival)  of  the  direct  wave  changes.  The  plots  are  given  for  each  of  the 
six  configurations  of  spheres  sketched  in  Figures  3-3,  3-5,  3-7,  3-9, 

3-11  and  3-13.  The  configurations  were  chosen  to  examine  the  varia¬ 
tions  as  (1)  the  number  of  spheres,  and  (2)  the  average  distance  of  the 
spheres  from  the  observation  point,  were  varied.  The  plots  for  con¬ 
stant  frequency  were  carried  out  at  the  geometric  mean  value  of  ft 
over  the  band  (rr/2  -  2  -t)  over  which  the  plots  for  constant  angle  of 
arrival  were  made.  The  plots  for  constant  angle  of  arrival  were 
carried  out  for  1  =  0. 

2 

The  average  value  of  the  computed  field  and  the  variance  a  of  all  the 

observations  are  computed  for  each  case  and  presented  in  the  figures. 

From  the  values  of  these  two  quantities,  it  appears  that  the  average 

value  of  the  amplitudes  for  changing  angle  and  frequency  are  more 

nearly  unity,  and  the  variances  for  the  two  cases  are  more  nearly 

equal,  for  a  large  number  of  spheres  at  four  or  more  wavelengths' 

average  distance  from  the  observation  point  (Figures  3-12  and  3-14). 

2 

Note  particularly  the  inverse -square  dependence  of  a  on  distance  for 

2 

eight  spheres;  this  is  what  the  statistical  analysis  predicts,  since  o 
would  then  be  half  the  average  scattered  power  at  the  observation  point. 

The  method  emplpyed  to  determine  whether  the  case  of  many  spheres, 
four  or  more  wavelengths  distant,  can  be  interpreted  statistically  is 
described  in  detail  in  Appendix  B.  By  means  of  this  method  the  proba¬ 
bility  that  the  fluctuations  in  amplitude  are  described  by  the  hypothesized 
distribution  may  be  calculated.  In  the  first  column  of  Table  3-1 
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FIGURE  3-7  CONFIGURATION  OF  FOUR  SPHERES  ARRANGED  ABOUT 
OBSERVATION  POINT  (LARGEST  MEAN  DISTANCE) 


FIGURE  3-12  TOTAL  (DIRECT  PLUS  SCATTERED)  ELECTRIC  FIELD  AMPLITUDE  (EIGHT  SPHERES) 


TOTAL  (DIRECT  PLUS  SCATTERED)  ELECTRIC  FIELD  AMPLITUDE  (  EIGHT  SPHERES) 


k2(angle)  \2 

probability 

k2(frequency)  y2 

probability 

.  0230 

>30% 

.  0180 

>50% 

.0256 

>20% 

.  0232 

>20% 

.  0198 

>10% 

.  0182 

>50% 

.0242 

>30% 

.  0184 

>50% 

.  0154 

>  5% 

.  0236 

>30% 

.  0204 

>20% 

.  0248 

>50% 

.  0262 

>30% 

.  0242 

>50% 

.  0236 

>50% 

.  0268 

>50% 

.  0244 

>  5% 

.  0588 

>10% 

.0210  Increasing 

>30% 

. 0172  Increasing 

>10% 

.0216  Frequency 

>  5% 

.  02  30  Angle 

>50% 

. 0244  1 

>50% 

.0374  l 

>50% 

.  0192 

>20% 

.  0378 

>70% 

.  0240 

>10% 

.  0284 

>50% 

.  0210 

>70% 

.  0280 

>50% 

.  0252 

>50% 

.  0268 

>20% 

.0256 

>20% 

.  0320 

>70% 

.  0204 

>  5% 

.  0298 

>50% 

.0262 

>20% 

.  0332 

>50% 

.02  78 

>20% 

.  0390 

>  5% 

.  0212 

>50% 

.  0360 

>50% 

.0324 

>  5% 

.  0328 

>50% 

.  0342 

>50% 

.  0294 

>30% 

.0340 

>10% 

.0176 

>50% 

.  0406 

>50% 

.  0332 

>50% 

.  0218 

>50% 

.  0290 

>20% 

.  0246 

>50% 

.  0288 

>70% 

.  0222 

>30% 

.  0254 

>50% 

.  0372 

>50% 

.  0204 

>50% 

.  0268 

>70% 

.  0272 

>50% 

.  0376 

>10% 

.  0244 

>30% 

.0322  . 

>50% 

.  0172 

>10% 

.0290 

>70% 

.  0142 

>70% 

.  0308 

>30% 

.  0224 

>70% 

.0328 

>50% 

.  0156 

>50% 

.  0306 

>50% 

.  0210 

>50% 

2 

Avg.  k  (angle)  = 

.  0262 

2 

Avg.  k  (frequency) 

=  . 0258 

2 

Std.  dev.  k  (angle)  = 

.  0058 

2 

Std.  dev.  k  (frequency)  =  .  0076 

TABLE  3-1 

k2  (ANGLE)  AND  k2  (FREQUENCY)  FOR 
36  VALUES  OF  FREQUENCY  AND  ANGLE  RESPECTIVELY 


are  given  36  values  of  k  (angle)  calculated  at  36  equis paced  values  of 
ft  from  ft  =  ~  to  ft  =  TT.  In  the  second  column  are  given  the  probabili¬ 
ties  that  the  independent  samples  over  angle  have  the  hypothesized 

2 

distribution,  as  measured  by  y  .  In  the  third  column  are  given  36 
2 

values  of  k  (frequency),  measured  at  36  equispaced  values  of  i|i  from 

l|l  =  0  to  t|r  =  2tt.  In  the  fourth  column  are  given  the  probabilities  that 

the  independent  samples  over  frequency  have  the  hypothesized  distri- 

2 

bution,  as  measured  by  y  .  The  calculations  pertaining  to  the 
configuration  of  eight  spheres  given  in  Figure  3-14  are  presented  in  a 
similar  fashion  in  Table  3-2. 


Several  comments  can  be  made  about  these  tables.  First  of  all,  as 

2 

indicated  by  the  values  of  y  probability,  the  data  are  in  good  agree¬ 
ment  with  the  hypothesized  distribution.  The  agreement  is  uniformly 

2 

good  for  the  various  values  of  k  (frequency),  and  appears  to  improve 

2 

somewhat  with  increasing  frequency  for  the  various  values  of  k  (angle) 

for  the  smaller  average  distance.  Secondly,  that  the  agreement  be- 

2 

tween  the  average  values  of  the  two  k  should  be  so  close  for  each 

table  is  not  too  surprising,  but  that  the  variation  in  the  values  of  the 
2 

two  k  should  be  comparatively  small  certainly  is.  The  largest  and 

2 

smallest  values  of  k  (angle)  in  Table  3-1  are  .0406  and  0154,  and  of 
2 

k  (frequency)  are  .0588  and  0142.  The  largest  and  smallest  values 

of  k^(angle)  in  Table  3-2  are  .  0104  and  .0030  and  of  k^(frequency)  are 

0116  and  .  0032.  This  indicates  that  a  single  measurement  at  one 

value  of  frequency  or  angle  is  likely  to  be  less  than  a  factor  of  two 

2  2 

incorrect  as  an  estimate  of  k  (angle)  or  k  (frequency)  respectively 

(see  the  discussion  given  in  Appendix  C).  Thirdly,  the  average  values 
2 

of  k  exhibit  a  strong  inverse  square  dependence  on  distance,  as  pre¬ 
dicted.  The  difference  between  the  computations  of  Table  3-1  and 


k2  (angle) 

probability 

k2  (frequency)  probability 

.0050 

>50% 

.  0042 

>70% 

.  0070 

>20% 

.  0068 

>50% 

.  0054 

>70% 

.  0052 

>70% 

.0030 

>30% 

.  0070 

>30% 

.0046 

>50% 

.  0050 

>50% 

.  0062 

>50% 

.  0038 

>70% 

.  0042 

>10% 

.  0064 

>70% 

.  0044 

>70% 

.  0084 

>30% 

.0038 

Increasing 

>30% 

.  0072 

Increasing  >70% 

.0070 

F  requency 

>70% 

.  0690 

Angle  >50% 

.  0050 

1 

>50% 

.  0106 

l  >  50% 

.0048 

>50% 

.  0052 

>90% 

.  0060 

>50% 

.  0066 

>70% 

.  0068 

>70% 

.  0016 

>90% 

.  0056 

>30% 

.  0082 

>50% 

.0074 

>50% 

.  0066 

>50% 

.  0070 

>50% 

.  0040 

>30% 

.  0044 

>70% 

.  0042 

>50% 

.  0076 

>30% 

.  0056 

>50% 

.  0068 

>30% 

.  0090 

>70% 

.0052 

>20% 

.  0082 

>70% 

.  0054 

>70% 

.  0060 

>70% 

.0060 

>50% 

.  0056 

>30% 

.0082 

>50% 

.  0058 

>70% 

.0092 

>90% 

.  0050 

>50% 

.  0084 

>10% 

.  0102 

>90% 

.0040 

>50% 

.  0074 

>70% 

.  0066 

>70% 

.  0072 

>70% 

.  0050 

>90% 

.  0064 

>50% 

.  0082 

>50% 

.  0054 

>50% 

.0062 

>50% 

.  0058 

>70% 

.0104 

>70% 

.  0044 

>90% 

.  0056 

>50% 

.  0032 

>30% 

.  0088 

>90% 

.  0044 

>70% 

.  0086 

>50% 

.  0052 

>70% 

.  0098 

>50% 

.  0062 

>50% 

Avg.  k2 

(angle)  = 

.  0063 

Avg.  k2 

(frequency)  =  .  0064 

Std.  dev.  k 

2  (angle)  = 

.0018 

Std.  dev.  k 

2 

(frequency)  =  .  0019 

TABLE  3-2 

k2  (ANGLE)  AND  k2  (FREQUENCY)  FOR 
36  VALUES  OF  FREQUENCY  AND  ANGLE  RESPECTIVELY 
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Table  3 -2  is  that  the  distance  from  each  of  the  spheres  to  the  observa- 

2 

tion  point  has  been  doubled,  and  the  average  values  of  k  are  observed 
to  drop  by  a  factor  of  four 

2  2 

Since  the  agreement  between  k  (angle)  and  k  (frequency)  is  good  for 

a  frequency  spread  of  two  octaves  it  seems  reasonable  to  inquire 

what  bandwidth  is  necessary  to  obtain  good  agreement.  To  this  end, 

2 

the  bandwidth  was  reduced  to  one  octave  and  k  (frequency)  was  meas¬ 
ured  in  each  of  the  two  cases  given  above.  Table  3-3  presents  the 
2 

value  of  k  (frequency)  for  18  equispaced  values  of  \|t  for  each  case. 

2 

The  average  values  of  k  (frequency)  arc  seen  to  be  in  good  agreement 

with  those  of  Tables  3-1  and  3-2.  Reducing  the  bandwidth  necessary 
2 

to  calculate  k  (frequency)  even  further  did  not  appear  practical,  since 
the  number  of  independent  sample  points  available  rapidly  diminishes 


The  resonance  region  was  chosen  as  the  frequency  range  of  most  in¬ 
terest,  but  the  statistical  model  proved  so  successful  in  this  range 
that  it  appeared  worthwhile  to  investigate  the  scattering  behaviour  of 
the  spheres  above  and  below  this  region  Accordingly,  brief  calcu¬ 
lations  were  made  of  the  relationship  between  k^(angle)  and  ^(fre¬ 
quency)  for  eight  spheres  for  val  ves  of  ft<TT/2  and  >2n.  It  was  found 

that  the  two  k^  were  approximately  equal  for  ^  in/4,  but  k^(angle) 

2 

started  to  become  significantly  less  than  k  (frequency)  below  this 

2 

value  of  h  It  was  found  that  k  (angle)  began  to  grow  appreciably 
larger  than  k^(frequency)  for  values  of  ft  2  3tt.  In  the  case  of  fi  2  rr/4 
it  appears  that  the  frequency  is  too  low  for  large  fluctuations  in  am¬ 
plitude  to  occur  as  the  angle -of-ar rival  is  varied  In  the  case  ol 
*  2  3tt,  the  large  forward -scattering  behaviour  of  the  spheres  caused 
great  changes  in  amplitude  to  occur  whenever  a  sphere  was  directly 
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ght  Spheres  Near  Observation  Point  Eight  Spheres  Far  From  Observation  Point 

2  2 
k  (frequency)  one  octave  k  (frequency)  one  octave 


0278 
0156 
.  0212 
.  0184 
0120 
0372 
0454 
.  0318 
0304 
0250 
0282 
0256 
.  0364 
0290 
0244 
0258 
0224 
0212 


0060 

0052 

0072 

0044 

0100 

0056 

0094 

0058 

0068 

0048 

0104 

0080 

0060 

0066 

0048 

0056 

0088 

0054 


Increasing 

Angle 

i 


Avg  k  (frequency)  =  0265 

2 

Std  dev  k  (frequency)  =  .  0074 


Avg  k  (frequency)  =  0067 

2 

Std  dev  k  (frequency)  =  .  0018 


TABLE  3-3 

k2  (ANGLE)  AND  k2  (FREQUENCY)  FOR 


18  VALUES  OF  FREQUENCY  AND  ANGLE  RESPECTIVELY 


between  the  direct  wave  source  and  the  observation  point  (i)r.  =  i|r  -  tt) . 

2  1 

Unless  the  cut  in  angle  taken  to  determine  k  (frequency)  also  coincides 
with  one  of  these  angles,  the  fluctuations  in  amplitude  with  frequency 
are  small  compared  to  the  large  forward-scattering  component. 

2  2 

Finally,  both  k  (angle)  and  k  (frequency)  were  computed  for  each  of 

10  configurations  of  spheres.  The  values  are  given  in  Table  3-4. 

Since  the  preceding  computation  has  indicated  that  a  one-octave  average 

2 

is  sufficient  to  provide  a  good  estimate  of  k  (frequency),  the  computa¬ 
tions  were  performed  for  one  octave  only.  The  values  given  in  the 

2 

table  are  averages  of  18  samples  of  k  (angle)  and  18  samples  of 
2 

k  (frequency).  Configurations  1  through  6  are  various  arrangements 
of  eight  spheres,  all  of  which  differ  from  the  two  previously  examined. 
In  configurations  1,  3,  and  5  the  closest  sphere  is  four  wavelengths 
distant  from  the  observation  point  at  the  lowest  frequency  at  which 
computations  were  made.  In  configurations  2,  4,  and  6  the  closest 
sphere  is  eight  wavelengths  distant  from  the  observation  point. 
Configurations  7  through  10  are  arrangements  of  16  spheres.  In  con¬ 
figurations  7  and  9  the  closest  sphere  is  four  wavelengths  distant.  In 

configurations  8  and  10  the  closest  sphere  is  eight  wavelengths  distant. 

2  2 

The  values  of  k  (angle)  and  k  (frequency)  for  each  of  the  eight-sphere 
configurations  are  higher  than  the  configurations  previously  considered, 
because  in  the  previous  configurations  the  larger  spheres  were  placed 
further  away  from  the  observation  point  than  the  smaller  spheres. 

The  values  given  in  Table  3-4  strongly  bear  out  the  hypothesis  that 

2  2 
the  two  k  are  equal  for  arrays  of  many  objects.  That  k  represents 

the  mean  scattered  power  at  the  observation  point  is  also  strongly 
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Configuration  No 

k  (angle) 

2 

k  (frequency) 

1 

0283 

0288 

2 

0068 

.  0067 

3 

.  0271 

.  0278 

4 

0065 

0065 

5 

0294 

.  0292 

6 

.  0072 

.  0071 

7 

0537 

.  0541 

8 

.  0127 

.0130 

9 

0552 

0560 

10 

0135 

0137 

TABLE  3-4 

AVERAGE  k2  (ANGLE)  AND  k2  (FREQUENCY)  FOR 
EACH  OF  10  CONFIGURATIONS  OF  SPHERES 
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borne  out  by  the  table,  as  for  each  configuration  k  is  seen  to  decline 

by  approximately  a  factor  of  four  when  the  mean  distance  from  the 

spheres  to  the  observation  point  is  doubled.  Another  indication  that 

2 

this  is  the  case  is  the  fact  that  the  value  of  k  approximately  doubles 
when  the  number  of  spheres  is  doubled. 

Calculation  of  electromagnetic  scattering  from  arrays  of  many 

spheres  have  therefore  shown  that  a  statistical  model  can  success  - 

2 

fully  represent  the  observed  effects.  In  addition,  the  values  of  k 
2 

(angle)  and  k  (frequency)  are  equal  for  a  large  number  of  spheres. 

2  2 
The  value  of  k  (frequency)  at  one  angle  or  of  k  (angle)  at  one  fre¬ 
quency  is  a  good  (within  a  factor  of  two)  estimate  of  the  many-angle 

average  or  many-frequency  average  oi  these  quantities.  Finally,  for 

2 

the  configurations  considered,  a  one -octave  measurement  of  k  (fre¬ 
quency)  is  sufficient. 

C.  Summary 

This  chapter  enlarges  the  scope  of  the  analytical  model  of  a  D/F  site 
by  suggesting  and  then  investigating  a  new  concept,  namely,  that  the 
characteristics  of  a  site  might  be  determined  by  averaging  over  fre¬ 
quency  rather  than  angle -of-arrival.  In  Section  A  the  reasonableness 
of  the  statistical  assumptions  applied  to  this  case  is  discussed.  In 
Section  B  the  concept  is  found  to  be  true  for  arrays  of  hemispherical 
bosses  on  a  conducting  ground  plane. 
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Chapter  4 


EXPERIMENTAL  PROCEDURE 


A.  Experimental  Model 

Before  describing  the  experimental  procedure  in  detail,  a  brief  intro¬ 
duction  is  given:  The  purpose  of  the  measurement  is,  basically,  to 
answer  two  questions  for  cases  which  cannot  be  done  theoretically. 

Does  the  analytical  model  fit  a  variety  of  sites,  both  for  varying  angle 

2 

of  arrival  and  for  varying  frequency?  Under  what  conditions  are  k 
2 

(angle)  and  k  (frequency)  approximately  equal?  In  order  to  answer 
these  questions,  some  means  must  be  found  to  measure  the  amount  of 
fluctuation  of  the  electric  field  magnitude  at  the  center  of  an  obstacle- 
strewn  D/F  site  as  (1)  a  plane  wave  of  constant  magnitude  and  frequen¬ 
cy  impinges  on  the  site  at  various  angles  of  arrival  or  (2)  a  plane  wave 
of  constant  magnitude  and  angle  of  arrival  but  variable  frequency 
impinges  on  the  site.  The  experimental  procedure  followed  in  the  first 
case  is,  in  essence,  to  move  a  target  transmitter  of  constant  frequen¬ 
cy  about  an  omnidirectional  (in  azimuth)  receiving  antenna  on  a  circle 
of  constant  radius  centered  on  the  receiving  antenna.  In  the  second 
case,  the  target  transmitter  is  fixed  in  azimuth  and  the  frequency  of 
transmission  is  swept  over  some  arbitrary  band.  In  each  case,  the 
observed  fluctuations  in  received  signal  level  are  proportional  to  the 
fluctuations  in  the  electric  field  magnitude  at  the  receiving  antenna. 
These  measured  fluctuations  are  then  analyzed  in  a  specific  way  to 
answer  the  questions  posed  earlier. 
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Checking  the  theoretical  conclusions  in  the  1  -  30  MHz  range,  with  a 
large  variety  of  sites,  presents  many  practical  difficulties.  Not  only 
is  a  sufficiently  large  variety  of  sites  unavailable  in  the  Syracuse  area, 
but  obstacles  of  the  requisite  size  and  number  are  not  readily  trans¬ 
portable.  In  addition,  the  variation  in  received  signal  level  with  angle 
of  arrival  is  not  easily  investigated,  if  the  radius  of  the  circle  along 
which  the  target  transmitter  must  move  is  large.  The  problems 
involved  in  keeping  to  a  precise  circular  path  two  miles,  say,  in 
diameter  are  obvious. 

In  order  to  simplify  the  experiment  as  much  as  possible  the  work  was 
carried  out  at  microwave  frequencies.  The  same  type  of  scattering 
phenomena  observable  in  the  HF  band  will,  with  a  suitable  reduction 
in  the  size  of  the  scatterers,  appear  in  the  microwave  region.  Also, 
if  the  electrical  parameters  of  the  media  used  in  the  experiment  are 
judiciously  chosen  to  match  those  of  the  full  size  media  in  the  HF 
band,  the  magnitude  of  the  scattered  waves  relative  to  the  directly 
transmitted  wave  is  the  same  order  of  magnitude  in  both  cases.  Spe¬ 
cifically,  if  the  complex  relative  dielectric  constant  of  the  material 
chosen  to  represent  the  earth  at  microwave  frequencies  is  in  the  range 
5  £  £  30,  3  u  e^'  £  30,  which  is  the  range  of  the  omplex  relative 

dielectric  constant  of  the  earth  at  HF,  then  the  material  is  suitable 
for  an  approximate  simulation  of  the  earth.  (It  should  be  noted  here 
that  this  statement  is  equivalent  to  saying  that  the  conductivity  of  the 
earth  by  a  factor  equal  to  the  ratio  of  operating  frequencies).  A  new 
type  of  carbon-loaded  dielectric  foam  manufactured  by  the  Emerson 
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and  Cuming  Company  undei  the  trade  name  of  Eccosorb  LS-26  has, 
at  microwave  frequencies,  a  complex  relative  dielectric  constant  of 
5  +  j  1 Z ,  so  it  meets  this  requirement. 

Figure  4-1  is  a  diagram  of  the  experimental  apparatus  and  Figure  4-2 
is  a  photograph  of  it.  It  consists  of  an  8'  x  12'  x  1/2"  plane  of  plywood 
on  which  is  placed  a  3/4"  thick  layer  of  LS-26.  Two  feet  from  one  end 
of  the  plywood  sheet  is  a  X/4  (at  5  GHz)  monopole,  used  as  the  tr  ins  - 
mitting  antenna.  The  monopole  is  constructed  of  a  brass  sleeve  fitted 
over  the  stub  end  of  a  type  N  panel  mount  (Figure  4-3).  At  the  other 
end  of  the  plywood  sheet  a  hole  4'  in  diameter  has  been  cut.  A  circular 
disk  of  3/4"  thick  pressboard  slightly  less  than  4'  in  diameter  is  cen¬ 
tered  in  the  hole  and  rotates  freely  within  it  on  eight  aluminum  rollers 
equispaced  around  the  periphery  of  the  disk.  At  the  center  of  the  press  - 
board  disk  is  another  X/4  monopole,  used  as  the  receiving  antenna. 

The  pressboard  is  also  covered  with  a  3/4"  thick  layer  of  LS-26.  The 
receiving  antenna  is  thus  placed  on  a  rotatable  "site"  approximately 
20\  (at  5  GHz)  in  diameter.  The  transmitting  antenna  is  fixed,  but 
rotating  the  site  around  the  receiving  antenna  is  equivalent  to  holding 
the  site  fixed  and  moving  the  transmitting  antenna  around  it,  as  long 
as  the  effect  of  only  those  objects  less  than  10  wavelengths  distant 
from  the  direction-finder  is  of  interest. 

During  the  measurement  procedure  various  obstacles  are  placed  on 
the  "site.  "  They  consist  of  either  various  pieces  of  LS-26  cut  to 
resemble  hills,  ridges,  etc.  ,  or  groups  of  thin  vertical  conductors  to 
represent  trees,  poles,  etc.,  or  metal  frameworks  to  represent 
buildings,  storage  facilities,  etc.  The  fluctuations  in  received  signal 


FIGURE  4-1  DIAGRAM  OF  EXPERIMENTAL  APPARATUS 
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FIGURE  4-2  PHOTOGRAPH  OF  EXPERIMENTAL  APPARATUS 
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level  are  observed  as  a  function  of  the  angular  position  of  the  press- 
board  disk  and  then  of  the  frequency  of  transmission  for  many 
different  configurations  of  these  obstacles. 

The  plywood  plane  rests  on  saw-horses  four  feet  above  the  floor  of 
the  anechoic  chamber  in  which  the  entire  assembly  is  placed.  The 
anechoic  chamber  functions  adequately  above  2  GHz.  The  pressboard 
disk  is  rotated  by  a  Scientific -Atlanta  azimuth  positioner,  which  is 
controlled  by  a  servomechanism  drive,  and  synchronized  with  a  polar 
recorder. 

A  block  diagram  of  the  transmitting- receiving  system  is  given  in 
Figure  4-4.  The  transmitting  system  consists  of  a  sweep  oscillator 
which  drives  a  microwave  amplifier.  The  oscillator  is  square-wave 
modulated  at  1  kHz.  The  amplifier  will  deliver  up  to  500  mW  of  power 
to  the  X/4  antenna.  The  transmitting  system  is  placed  beneath  the 
plywood  plane  to  minimize  interference  with  the  measurements. 

The  receiving  system  employs  a  funnel  diode  detector,  the  audio  out¬ 
put  of  which  is  fed  into  the  bolometer  amplifier  of  the  polar  recorder. 
The  tunnel  diode  detector  is  employed  instead  of  a  conventional  crystal 
detector  because  the  received  signal  level  is  on  the  order  of  -35  dBm. 
This  level  is  approximately  equal  to  the  tangential  sensitivity  of  a 
conventional  detector.  The  tangential  sensitivity  of  a  tunnel  diode 
detector,  on  the  other  hand,  is  approximately  -55  dBm. 

The  data  collection  system  functions  in  the  following  manner:  The  AC 
output  of  the  bolometer  amplifier  is  fed  into  a  high  input  impedance 
full-wave  rectifier.  The  DC  output  of  the  rectifier  is  sampled  by  a 
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FIGURE  4-4  BLOCK  DIAGRAM  OF  TRANSMITTING- RECEIVING  SYSTEM 


digital  voltmeter.  The  digitized  output  of  the  voltmeter  is  fed  into  a 
printer.  The  sampling  rate  of  the  voltmeter  is  synchronized  with  the 
rotation  rate  of  the  polar  recorder  to  provide  approximately  one 
sample  per  degree  of  angle.  When  the  variations  with  respect  to  fre¬ 
quency  are  investigated,  the  positioner  is  disconnected  from  the 
pressboard  disk,  and  the  rotation  rate  of  the  polar  recorder  is  syn¬ 
chronized  with  the  sweep  rate  of  the  oscillation.  In  this  case,  the 
angular  variable,  8,  is  directly  proportional  to  frequency.  The 
sampling  rate  of  the  voltmeter  is  then  adjusted  to  provide  approxi¬ 
mately  one  sample  every  10  MHz  of  frequency. 

The  experimental  apparatus  supplies  outputs  in  each  of  two  forms. 

The  polar  recorder  graphs  (linearly)  received  power  versus  angle  and 
versus  frequency.  The  digital  sampling  system  in  effect  samples 
these  graphs  and  provides  a  printed  list  of  these  samples.  Since  only 
relative  magnitudes  are  of  importance,  no  absolute  calibration  of  the 
system  is  necessary.  The  plotted  graphs  are  of  interest  only  as  a 
pictorial  indication  of  the  amount  of  variation  observed.  Several  of 
these  are  exhibited  in  Section  C  of  this  chapter.  All  of  the  analysis 
is  done  on  the  digital  samples  which  are  transferred  to  punched  cards 
and  analyzed  by  computer. 

B.  Adjustment  of  Data 

A  preliminary  look  at  the  data  indicated  that  it  would  have  to  be  ad¬ 
justed  before  analysis.  It  was  found  that  the  magnitude  of  the  direct 
wave  was  not  independent  of  angle  or  frequency.  This  was  determined 
from  measurements  made  without  any  obstacles  present  on  the  rotating 
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disk.  The  variation  with  angle  was  first  thought  due  to  mechanical 
problems  with  the  rotatable  site,  but  when  these  were  overcome,  some 
variation  persisted.  This  variation  is  ascribed  to  inhomogeneity  in 
the  Eccosorb  LS-26. 

The  manufacturer  states  that  reasonable  tolerances,  within  a  2'  x  2' 

sheet  of  the  material,  are  e1  =  5  ±  2,  e"  =  12  ±  5.  Differences  of 

r  r 

this  order  of  magnitude  are  sufficient  to  explain  the  0.  5  -  1  dB  variation 
in  the  observed  magnitude  of  the  direct  wave.  The  variation  with  fre¬ 
quency  is  a  result  of  several  factors: 

1.  The  antennas  are  \/4  long  at  the  mean  frequency  of 
the  band.  The  sensitivity  of  the  transmitting-receiving 
system  therefore  varies  with  frequency. 

2.  The  Eccosorb  is  dispersive.  Its  dielectric  constant 
varies  slowly  with  frequency. 

3.  The  output  of  the  sweeper-amplifier  combination  varies 
slightly  over  the  band,  regardless  of  how  well  leveled 
the  sweeper  is. 

Fortunately,  however,  the  variation  is  slow  compared  to  the  rapid 

fluctuations  of  the  scattered  wave  magnitude  E  .  This  suggests  that 

s 

a  moving  average  be  computed  and  the  ratio  of  the  amplitude  of  the  total 
(direct  plus  scattered)  signal  to  the  amplitude  of  the  moving  average 
signal  be  used  as  the  variable  r.  (Note  that  this  procedure  is  necessary 
in  the  full-scale  case  as  well,  since  it  is  expected  that  the  amplitude 
of  the  direct  wave  will  be  similarly  dependent  on  the  path  attenuation 
and  the  type  of  transmitting-receiving  system  used.  )  The  moving 
average  is  calculated  in  order  to  remove  deterministic  quantities  from 
the  data.  (See  Appendix  C) 
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Figure  4-5  shows  a  moving  average  computed  over  45°  (dotted  line) 
and  the  actual  measurement  with  obstacles  present  (solid  line). 

Figure  4-6  shows  a  moving  average  computed  over  400  MHz  (dotted 
line)  and  the  measurements  with  obstacles  present  (solid  line).  The 
averaging  intervals,  45°  and  400  MHz,  respectively,  were  chosen  as 
those  which  most  accurately  reproduced  the  case  with  obstacles  absent 
(Figure  4-7  for  angle  and  Figure  4-8  for  frequency).  It  was  decided 
to  compute  this  moving  average  rather  than  remeasure  the  obstacle- 
free  case  each  time  because  this  is  what  one  would  do  in  an  actual, 
full-scale,  investigation  (the  obstacle-free  case  being  unavailable). 

In  the  figures,  the  ordinate  is  the  measured  amplitude  of  the  received 
voltage  normalized  to  its  average  value.  This  is  done  to  remove  the 
effect  of  gain  settings  on  the  presented  amplitudes. 

The  independent  sample  points  were  determined  according  to  the 
procedure  given  in  Appendix  B.  The  correlation  between  pairs  of 
points  was  calculated  and  the  interval  between  samples  was  chosen  as 
that  which  gave  zero  correlation.  The  independent  sampling  points 
were  then  analyzed  according  to  the  modified  minimum  method. 

C.  Experimental  Results 

The  obstacles  introduced  onto  the  site  were  of  several  types: 

1.  Pieces  of  Eccosorb  LS-26  cut  to  resemble  hills,  ridges,  etc. 

2.  Metal  frameworks  in  the  shape  of  buildings  or  storage 
facilities, 

3.  Clusters  of  vertical  conductors  to  represent  groups  of  trees. 


4.  Long  horizontal  conductors  to  represent  fences,  telephone 
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Measured  Amplitude  (Normalized) 
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FIGURE  4-7  MEASURED  AMPLITUDE  FOR  VARYING  ANGLE  CASE  -  NO  OBSTACLES 


It  should  be  emphasized  that  it  is  not  important  that  the  correspondence 
between  the  microwave  obstacles,  and  what  they  are  chosen  to  repre¬ 
sent,  be  exact.  The  microwave  experiment  stands  by  itself;  the 
analysis  given  in  Chapter  2  may  be  applied  to  it  regardless  of  whether 
or  not  it  is  a  precise  representation  of  the  HF  problem.  The  micro- 
wave  obstacles,  however,  have  been  chosen  to  be  similar  to  naturally 
occurring  obstacles  in  the  HF  band  in  order  to  support  the  application 
of  the  statistical  theory  to  HF.  It  is  certainly  reasonable  to  propose 
that  if  the  statistical  analysis  is  successful  in  the  microwave  region 
it  will  also  be  successful  at  HF  when  applied  to  obstacles  of  similar 
dimensions  (with  respect  to  a  wavelength)  and  similar  electrical 
characteristics. 

The  measurements  for  the  varying  angle  case  were  performed  at 
5  GHz,  The  measurements  for  the  varying  frequency  case  were  per¬ 
formed  over  an  octave  band  from  8  GHz  to  4  GHz.  This  frequency 
range  was  initially  chosen  because  it  is  the  range  over  which  the  sweep¬ 
er  operates.  Other  measurements  were  performed  with  an  increased 

range  (down  to  2  GHz  and  up  to  12  GHz),  but,  for  the  obstacle  con- 

2 

figurations  considered,  the  value  of  k  (frequency)  increased  very 
slowly  as  the  bandwidth  was  increased. 

Twenty-four  different  configurations  were  examined.  Figures  4-9, 

4-10,  4-11  and  4-12  are  the  polar  plots  of  received  signal  power 
versus  angle  of  arrival  for  each  configuration.  Figures  4-13,  4-14, 
4-15  and  4-16  are  the  polar  plots  of  received  signal  power  versus 
frequency  for  each  configuration.  There  are  six  configurations  per 
figure.  Each  plot  is  linear,  with  the  radius  proportional  to 
received  power. 
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Table  4-1  summarizes  the  results  of  the  measurements,  after  the  ad¬ 
justments  described  in  Chapter  2-B  are  performed.  Column  1  gives 

the  number  of  the  obstacle  configuration  presented.  Column  2  is  the 
2 

value  of  k  for  the  case  of  varying  angle.  Column  3  states  the 

2 

probability,  as  measured  by  y  •  that  the  data  has  the  hypothesized 

2 

distribution,  for  the  varying  angle  case.  Column  4  is  the  value  of  k 

for  the  case  of  varying  frequency.  Column  5  is  the  probability,  as 
2 

measured  by  y  »  that  the  data  has  the  hypothesized  distribution,  for 

2 

the  varying  frequency  case.  Column  6  is  the  ratio  of  k  for  the 

2 

varying  angle  case  to  k  for  the  varying  frequency  case.  Column  7  is 
2 

the  value  of  k  (angle)  for  various  cases  normalized  to  the  value  of 
2 

k  (angle)  for  the  first  case  measured  with  obstacles  present.  Column 
2 

8  is  the  value  of  k  (frequency)  for  various  cases  normalized  to  the 
2 

value  of  k  (frequency)  for  the  first  case  measured  with  obstacles 
present. 


A  description  of  each  obstacle  configuration,  which  appears  in  the 

plots,  and  the  table,  follows.  A  few  representative  sketches  are 

given  in  Figures  4-17  and  4-18.  For  the  sake  of  convenience,  they 

are  referred  to  according  to  what  they  are  meant  to  represent  at  HF 

instead  of  what  they  are  (e.  g.  ,  "hill"  instead  of  "piece  of  Eccosorb"). 

The  configurations  were  chosen  so  that  the  effect  of  three  quantities 
2 

on  the  value  of  k  could  be  determined:  1)  the  number  of  obstacles 
present;  2)  the  distance  from  the  direction-finder  to  the  obstacles; 

3)  the  type  of  obstacles.  In  the  listing,  the  wavelength  X  refers  to 
the  wavelength  at  5  GHz. 

1.  No  obstacles  present. 

2.  Two  irregularly  shaped  hills,  approximately  X /4  high  and 
X  / 2  -  X  long,  placed  >2X  away  from  the  receiving  antenna. 
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(ANGLE)  AND  k  (FREQUENCY)  FOR  VARIOUS  OBSTACLE  CONFIGURATIONS 


ROTATING  DISC 


FIGURE  4-17  OBSTACLE  ARRANGEMENT  INCASE 


ROTATING  DISC 


FIGURE  4-18  OBSTACLE  ARRANGEMENT  IN  CASE  16 


3.  Four  irregularly  shaped  hills,  of  the  same  general  character¬ 
istics,  placed  >2X  away  from  the  receiving  antenna. 

4.  Eight  irregularly  shaped  hills,  of  the  same  general  character¬ 
istics,  placed>2X  away  from  the  receiving  antenna. 

5.  Two  of  the  same  obstacles,  placed  >4X  away  from  the  re¬ 
ceiving  antenna. 

6.  Four  of  the  same  obstacles,  placed  >4X  away  from  the  re¬ 
ceiving  antenna. 

7.  Eight  of  the  same  obstacles,  placed  >4X  away  from  the  re¬ 
ceiving  antenna. 

8.  Two  of  the  same  obstacles,  placed  >6X  away  from  the  re¬ 
ceiving  antenna. 

9.  Four  of  the  same  obstacles,  placed  >6X  away  from  the  re¬ 
ceiving  antenna. 

10.  Eight  of  the  same  obstacles,  placed  >6X  away  from  the  re¬ 
ceiving  antenna. 

11.  A  clump  of  trees  <X/4  high  and  ~X  long,  and  a  building  <X/4 
high  and  ~X  long,  placed  >2X  from  the  receiving  antenna. 

12.  A  telephone  line  ~X  / 4  high,  and  a  hill  <X  /4  high  and  ~X 
long  added  to  the  obstacles  in  1 1  and  placed  >2X  from  the  re¬ 
ceiving  antenna. 

13.  Two  clumps  of  trees,  two  buildings,  three  hills,  and  a  tele¬ 
phone  line,  all  with  the  same  characteristics  as  in  11  and  12, 
and  all  placed  >2X  from  the  receiving  antenna. 
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14.  The  same  obstacles  as  in  11,  placed  >4\  from  the  receiving 
antenna. 

15.  The  same  obstacles  as  in  12,  placed  >4\  from  the  receiving 
antenna. 

16.  The  same  obstacles  as  in  13,  placed  >4\  from  the  receiving 
antenna. 

17.  The  same  obstacles  as  in  11,  placed  >6X  from  the  receiving 
antenna. 

18.  The  same  obstacles  as  in  12,  placed  >6X  from  the  receiving 
antenna. 

19.  The  same  obstacles  as  in  13,  placed  >6\  from  the  receiving 
antenna. 

20.  Twenty  obstacles,  of  various  sizes  and  shapes,  all  <\/4  high 
and  between  X  / 4  and  2X  long,  and  placed  >2\  from  the  re¬ 
ceiving  antenna. 

21.  The  same  20  obstacles,  placed  >4X  from  the  receiving 
antenna. 

22.  As  in  21,  the  obstacles  rearranged. 

23.  The  same  20  obstacles,  placed  >6\  from  the  receiving 
antenna. 

24.  As  in  23,  the  obstacles  rearranged. 

The  results  in  Table  4-1  call  for  some  comment.  A  list  of  observations 
on  the  data  is  given,  with  some  explanations  of  the  observed  effects. 
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Using  the  5  per  cent  criterion  for  y  »  all  the  configurations 

given  in  the  table  fall  within  the  hypothesized  distribution. 

Surprisingly,  even  situations  with  only  two  obstacles  have  a 

high  probability  of  being  in  this  class,  indicating  that  E  cos  $ 

n  n 

and  E  sin  I  (n  =  1,  2)  are  sometimes  close  to  a  normal  dis- 
n  n 

tribution.  Since  most  of  the  probabilities  given  in  the  table 
are  significantly  greater  than  the  minimum  necessary,  the 
analytical  model  appears  to  explain  successfully  the  measured 
fluctuations  in  received  power. 

2  2 

The  value  of  k  (angle)  and  k  (frequency)  increases  as  the 
number  of  obstacles  increases.  This  is  what  is  expected, 
since  each  additional  obstacle  is  another  source  of  scattered 
power. 

2 

The  value  of  k  ,  for  both  cases,  decreases  as  the  obstacle- 

free  radius  around  the  receiving  antenna  increases.  The  value 
2 

of  k  is  approximately  inversely  proportional  to  the  distance 

from  the  obstacles  to  tlfb  receiving  antenna.  The  relationship 
2 

between  k  and  this  distance  is  difficult  to  establish,  how¬ 
ever,  since  in  moving  the  obstacles  outward  from  the  receiving 
antenna,  their  relative  orientation  cannot  be  maintained  easily. 
Thus  another  variable,  whose  influence  is  not  known,  is 
unavoidably  introduced. 

2  2 

The  value  of  k  (angle)  is  greater  than  k  (frequency)  for  all 

2 

configurations  considered.  As  previously  noted,  k  (frequen- 

2 

cy)  appears  to  converge  very  slowly  to  k  (angle)  as  the 
bandwidth  of  the  sweep  increases. 
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5.  The  agreement  between  k  (angle)  and  k  (frequency)  im¬ 
proves  as  the  number  of  obstacles  increases,  becoming  very 
good  for  a  large  number  of  obstacles.  This  may  result  from 

the  fact  that  the  number  of  independent  samples  of  the  data, 

2 

and,  correspondingly,  the  accuracy  in  the  estimate  of  k  , 
increases,  as  the  number  of  obstacles  increases,  for  both 
the  varying  angle  and  the  varying  frequency  case.  The 
number  of  independent  samples  is  roughly  dependent  on  the 
rapidity  of  the  fluctuations  in  the  received  power.  These 
fluctuations  are  much  more  rapid  as  the  number  of  obstacles 
increases. 

.  2  2 

6.  The  agreement  between  k  (angle)  and  k  (frequency)  im¬ 
proves  as  the  obstacle-free  radius  around  the  receiving 
antenna  increases.  This  is  traceable  to  the  same  sources  as 
5,  above,  since  the  rapid  "y  of  the  fluctuations  in  the  received 
power  increases  as  the  distance  to  the  obstacles  increases 
and  hence  more  independent  samples  of  the  data  are  available. 

2  2 

From  these  observations,  either  k  (angle)  or  k  (frequency)  appear 

to  satisfy  the  requirements  that  one  would  intuitively  expect  a  "site 

2 

goodness"  parameter  to  meet.  If  increasing  k  is  taken  to  mean  that 
a  site  is  less  suitable  for  direction-finding,  then  the  expected  cases 
exhibit  this  behavior,  namely: 

1.  The  closer  the  obstacles  are,  the  more  inaccuracy  is  intro- 

2 

duced  into  the  D/F  measurement,  k  increases  as  the 
obstacles  are  moved  closer. 
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2.  The  more  obstacles  there  are,  the  more  the  accuracy  of  the 

2 

D/F  system  is  impaired,  k  increases  as  the  number  of 
obstacles  increases. 

3.  The  larger  the  obstacles  are,  the  more  error  they  introduce. 

increases  as  the  obstacle  size  increases. 


2  2 
Since  k  (frequency)  exhibits  the  same  trends  as  k  (angle),  as  columns 

2 

7  and  8  of  Table  4-1  show,  it  is  just  as  useful  a  parameter  as  k  (angle) 

2 

even  though  there  is  a  discrepancy  between  it  and  k  (angle)  for  finite 
bandwidths.  In  fact,  it  is  probably  unnecessary  to  increase  the  band¬ 
width  in  practical  situations,  since  various  sites  may  still  be  compared 

2 

with  each  other  on  the  basis  of  a  "standard  bandwidth"  k  (frequency). 
The  following  section  will  elaborate  on  how  site  comparisons  on  this 
basis  might  be  accomplished. 


D.  Summary 


The  experimental  measurements  carried  out  to  verify  the  validity  of 

the  analytical  model  are  presented  in  this  Chapter.  In  Section  A  is 

given  a  description  of  the  experimental  apparatus,  the  transmitting- 

receiving  system,  and  the  automated  data  collection  system.  The 

experiments  were  carried  out  at  microwave  frequencies  for  reasons 

of  convenience,  but  because  of  the  electrical  parameters  chosen,  the 

observed  scattering  effects  are  similar  to  those  expected  in  the  HF 

band.  A  necessary  adjustment  of  the  data  is  presented  in  Section  B. 

2 

In  Section  C  the  experimental  results  are  given.  The  value  of  k 
2 

(angle)  and  k  (frequency)  is  calculated  for  each  configuration  of  obsta¬ 
cles  tested  and  the  probability  that  such  a  configuration  is  described 


by  the  analytical  model  is  given.  In  addition  several  observations  are 

2  2 

made  concerning  the  nature  of  k  {angle)  and  k  (frequency)  and  their 
usefulness  as  measures  of  site  suitability. 


Chapter  5 

CONCLUSIONS  AND  RECOMMENDATIONS 

A.  Conclusions 

The  commonly  accepted  method  of  D/F  site  evaluation*  is  to  move  a 
target  transmitter  of  fixed  frequency  on  a  circular  path  of  large  (with 
respect  to  a  wavelength)  radius  centered  on  the  direction-finder.  If 
the  maximum  variation  of  the  received  signal  as  a  function  of  the  an¬ 
gular  position  of  the  target  transmitter  is  greater  than  some  arbitrary 

•  •  ,v 

criterion,  the  site  is  rejected  for  direction-finding.  Since  it  is  usually 
desired  to  use  the  direction-finder  over  a  frequency  band,  this  test  is 
repeated  for  a  number  of  frequencies  within  the  band. 

It  is  possible  to  represent  in  a  three-dimensional  plot  the  electric 
field  amplitude  at  the  direction-finder  as  a  function  of  the  angular 
position  and  the  frequency  of  the  target  transmitter.  The  amplitude  is 
the  z-coordinate  of  a  surface  whose  x  and  y  coordinates  are  angle  and 
frequency  respectively.  The  method  described  above  gives  the  value 
of  z  along  several  lines  of  constant  y.  The  surface  could  therefore 
be  completely  described  if  a  sufficient  number  of  these  lines  were 
chosen.  Similarly  the  surface  could  be  completely  described  if  the 
value  of  z  were  given  along  a  sufficient  number  of  lines  of  constant  x. 

In  other  words,  the  evaluation  method  described  in  the  first  paragraph 
can  be  carried  out  equivalently  by  placing  the  target  transmitter  at  a 
number  of  points  on  the  chosen  circular  path  and  varying  the  frequency 
over  the  band  of  interest,  a  seemingly  trivial  observation.  In  practice, 
however,  this  latter  measurement  is  easier  to  perform. 
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The  difficulty  involved  in  the  former  measurement  lies  in  keeping  the 
radius  of  the  circle  constant  for  large  radii.  This  is  necessary  to 
insure  that  the  observed  perturbations  in  the  received  signal  arise  not 
from  variations  in  the  transmitter-receiver  distance,  but  only  from 
the  fluctuations  in  the  amplitudes  and  phases  of  the  scattered  waves 
brought  about  by  the  changing  angle  of  arrival  of  the  direct  wave.  In 
the  latter  case,  the  changes  in  the  absolute  magnitude  of  the  direct 
wave  caused  by  variations  in  the  transmitter-receiver  distance  do  not 
affect  the  result.  The  quantity  of  interest  is  the  maximum  fluctuation 
from  the  average  signal  magnitude  over  the  band,  where  it  is  assumed 
that  allowance  is  made  for  the  slowly-varying -with-frequency-sensi- 
tivity  of  the  transmitting -receiving  system. 

The  above  observation  is  true  of  any  site.  The  question  of  most  in¬ 
terest  to  the  engineer,  though,  is  how  many  points  on  the  approximate 
circular  path  must  be  used  to  evaluate  the  site.  There  is  no  single 
answer  to  this  question  applicable  to  all  sites.  In  the  body  of  this  work, 
however,  it  is  shown  that  for  a  certain  class  of  sites  very  few,  perhaps 
only  one  or  two,  points  are  necessary.  In  general  these  sites  are 
characterized  by  the  presence  of  several  "resonant"  obstacles;  i.  e.  , 
obstacles  whose  dimensions  are  on  the  order  of  a  wavelength  in  magni¬ 
tude  in  the  frequency  band  of  interest. 

For  this  type  of  site  a  one -parameter  probability  distribution  can  be 
derived  for  the  electric  field  amplitude  at  the  direction  finder  as  either 
(A)  the  angle  of  arrival  or  (B)  the  frequency  of  the  direct  wave  is 
varied.  The  computations  described  in  Chapter  3  and  the  experimen¬ 
tal  work  given  in  Chapter  4  indicate  that  the  parameter  of  case  (A)  is 
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approximately  the  same  as  the  parameter  of  case  (B)  when  several 

resonant  obstacles  are  present.  In  addition,  the  results  of  Chapters  3 

and  4  strongly  suggest  that  only  a  small  number  of  target  transmitter 

positions  are  necessary  to  provide  sufficient  data  to  calculate  this 

parameter  when  the  technique  of  frequency  variation  is  employed.  This 
2 

parameter,  k  ,  which  is  the  ratio  of  the  mean  square  power  in  the 
scattered  waves  to  the  power  in  the  direct  wave,  is  proposed  as  a 
figure  of  merit  for  sites  having  the  characteristic  mentioned  above. 

The  computational  work  has  also  given  indications  of  when  a  small 

number  of  transmitter  positions  may  not  be  sufficient  to  evaluate  a 

site.  Above  the  resonance  region,  for  example,  extended  obstacles 

are  more  directional  reradiators  than  in  the  resonance  region.  An 
2 

estimate  of  k  by  frequency  variation  in  this  range  is  therefore  criti¬ 
cally  dependent  on  the  transmitter's  angular  position.  Below  the 
resonance  region  the  magnitude  of  the  fluctuations  in  the  scattered 

power  as  the  angular  position  of  the  transmitter  is  changed  is  over- 
2 

estimated  by  a  k  measured  through  the  resonance  region. 

The  question  of  how  many  obstacles  must  be  present  for  a  rapid  deter- 
2 

mination  of  k  by  frequency  variation  has  been  left  unanswered.  The 
computational  and  experimental  work  has  indicated  that  the  answer  to 
this  question  depends  greatly  on  the  types  of  obstacles  involved;  i,  e.  , 
their  shape,  conductivity,  etc.  A  single  obstacle  with  a  complex  geo¬ 
metrical  form  may  satisfy  the  requirements  as  well  as  many  obstacles 
with  simple  geometrical  shapes. 
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Finally,  the  experimental  confirmation  is  somewhat  limited.  Although 

2 

in  the  computational  work  it  was  indicated  that  the  k  estimated  by 

2 

frequency  variation  would  be  different  from  the  k  estimated  by  angle  - 
of-arrival  variation  for  obstacles  large  with  respect  to  a  wavelength, 
an  extensive  investigation  of  this  effect  was  not  undertaken  in  the  ex¬ 
perimental  work.  Thus  it  is  difficult  to  define  precisely  the  obstacle 

2 

size  beyond  which  k  may  not  be  determined  by  a  small  number  of 

measurements.  For  the  range  of  obstacle  sizes  used  in  the  experimen- 

2 

tal  work,  however,  the  k  estimated  by  frequency  variation  at  a  single 
transmitter  position  is  very  close  to  that  determined  by  a  'gle -of- 
arrival  variation  at  the  geometric  mean  frequency  of  the  band. 


From  the  table  given  in  Chapter  4  it  can  be  seen  that  if  all  the  k  (angle) 

2  2 

are  normalized  to  one  value  of  k  (angle)  and  similarly  all  the  k  (fre¬ 
quency)  are  so  normalized  (columns  7  and  8),  then  each  normalized 

2  2 
k  (angle)  is  within  a  factor  of  two  of  each  normalized  k  (frequency) 

for  all  configurations  except  number  23.  Indications  are  strong 
2  2 

therefore  that  k  (frequency)  =  k  (angle)  except  for  a  bandwidth- 
dependent  multiplicative  factor. 


B.  Recommendations 

Several  recommendations  for  futher  work  in  this  area  have  evolved 
from  the  investigations  reported  here.  A  list  of  them  follows: 

1.  The  investigation  was  concerned  with  the  use  of  the  statistical 
methods  employed  in  the  analysis  of  rough- surface  scattering 
to  provide  a  rational  basis  of  D/F  site  selection.  No  attempt 
has  so  far  been  made  to  suggest  a  means  of  correcting  for  the 
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errors  introduced  by  the  scatterers.  The  previous  work 
mentioned  in  the  introduction  outlines  a  means  of  correcting 
for  the  presence  of  one  point  scatterer.  The  method  given  in 
Appendix  D  is  a  more  general  method  to  correct  for  the  errors 
introduced  by  several  point  scatterers,  but,  as  can  be  seen 
from  the  complexity  of  the  mathematics,  this  solution  becomes 
quite  cumbersome  for  more  than  three  scatterers.  In  addi¬ 
tion,  it  places  severe  requirements  on  the  antenna  pattern  of 
the  receiving  antenna  Accordingly,  it  is  recommended  that 
further  work  be  carried  out  to  provide  simpler  means  of 
correcting  for  the  presence  of  a  number  of  scatterers  during 
the  D/F  procedure. 

The  fluctuations  in  received  signal  level  versus  varying  ele- 
vational  angle  of  arrival  could  be  examined  theoretically  and 
experimentally.  The  assumptions  made  in  Chapter  2  may 
apply  to  this  case  as  well. 

Further  analytical  investigation  of  the  surface  and  volume 

current  distributions  on  and  in  a  complex  obstacle  is 

suggested.  Although  such  an  investigation  must  necessarily 

be  approximate,  it  would  serve  to  cast  further  light  on  what 

2 

conditions  must  be  met  in  order  that  k  (angle)  be  the  same 
2 

as  k  (frequency)  for  a  single  obstacle. 

Full  scale  experimental  work  is  suggested,  preferably  in 

a  geographic  area  where  a  good  site  is  readily  available. 

2 

The  change  in  the  value  of  k  as  complex  obstacles  are 
added  to  the  site  may  then  be  observed. 


5.  The  problem  of  siting  a  hydrophone  array  on  the  bottom  of  the 
ocean  is  analogous  to  that  of  siting  ^  radio  direction-finder. 
Random  scattering  of  sound  waves  from  various  solid  obstacles 
and  eddies  disrupt  the  performance  of  such  an  array.  The 
statistical  analysis  could  be  restated  in  acoustical  terms  and 
applied  to  this  problem  as  well. 

6.  The  angle -frequency  relationship  noted  in  this  work  may  have 
an  analogy  in  a  time-frequency  relationship  in  ionospheric 
propagation.  Perhaps  the  statistical  characteristics  of  the 
ionosphere  can  be  determined  by  frequency  variation  also. 

C.  Summary 

The  general  conclusions  arrived  at  from  the  study  and  some  recom¬ 
mendations  for  further  work  are  presented  in  this  chapter,  In  Section 
A  the  applicability  of  the  statistical  methods  is  discussed  and  a  simpli¬ 
fication  of  the  standard  method  of  site  evaluation  is  suggested.  In 
Section  B  further  work  in  this  area,  both  analytical  and  experimental 
is  recommended. 


Appendix  A 


Narrow  aperture  direction-finding  systems  typically  function  by  using 
figure  8  or  cardioid  antenna  patterns  to  directly  null  out  the  signal 
received  from  a  distant  transmitter.  The  depth  of  the  null  determines 
how  accurately  the  measurement  can  be  made,  as  illustrated  in  Fig¬ 
ure  A-l.  The  uncertainty  of  the  measurement  increases  as  the  depth 
decreases.  The  null  depth  is  affected  by  phase  and  amplitude  imbal¬ 
ances  in  the  antenna  structure,  a  low  signal-to-noise  ratio  of  the 
receiver,  polarization  error,  and  other  failures  of  the  system  to  per¬ 
form  according  to  its  theoretical  capacity.  It  is  also  affected  by 
external  noise  and  multipath  propagation  caused  by  site  nonuniform¬ 
ities  . 


The  parameter  k  is  the  average  power  contained  in  the  resultant  of 
the  scattered  waves  to  the  power  contained  in  the  direct  wave  and  may 
be  directly  related  to  the  depth  of  the  null.  For  examj#e,  for  two 
vertical  electric  dipoles  spaced  a  distance  d  apart,  and  connected 
180  out  of  phase  with  each  other,  the  voltage  response  in  the  azi¬ 
muthal  plane  to  the  direct  wave  is  given  by 

V  =  AE  sin  (  ~  sin  9  )  (Eq.  A-l) 

where  A  is  the  voltage  response  of  the  system  to  a  unit  field  and  9^ 
is  the  angle  of  arrival  of  the  direct  wave.  The  voltage  response  to 
the  scattered  waves  is  given  by: 


V 

s 


A 


N 

V 

f 

n=  I 


„  n  .  .  nd  . 

E  e  sin  (  — —  sm  9  ) 
n  >  n 


(Eq.  A-2) 
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where  8^  is  the  angle  of  arrival  of  the  nth  scattered  wave.  The  power 
received  from  the  scattered  waves  averaged  over  all  angles  of  arrival 
of  the  direct  wave  is: 


P  = 


<  V2> 
s 


£  £  j(*  -  *  ) 

:I  I  EnEme  "  m  sin  <  TT  Si"  V 
n  =  l  m=l 


sin  (  ~  sin  8  )>  (Eq.  A-3) 

A.  m 


where  Z  is  the  input  impedance  to  the  receiver.  Assume  E . E  , 

1  n 

i,  ....  $  are  independent  random  variables,  as  before,  and,  in  ad- 
1  n 

dition,  that  9  is  a  random  variable  independent  from  E  and  $  and 
n  r  n  n 

uniformly  distributed  from  0  to  2tt.  Then  (A-3)  becomes: 

„  N 


o  A  V  tr>2  .  2  ,  TTd  .  .  . 

P  =  ~  )  <E  >  <sin  (  —  sin  9  )> 
s  Z  L  n  X  n 

n  =  1 


(Eq.  A -4) 


Now 


<  sin 


2tt 

.  nd  .  .  1  ( 

(  —  sin  8  )>  =  —  \ 

X  n  2tt  j 


.  2  nd  .  «... 

sin  (  —  sin  9  )  d8 
X  n  n 


=  2  (*  ‘  Jo  (2?>)  forall9n  (Eq.A- 


5) 


so  we  can  write 


A2  2 

<P  >  =  “  <e2> 

s  Z  s 


(Eq.  A  -6) 
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If  the  power  contained  in  the  direct  wave  is  much  greater  than  the 

2 

power  contained  in  all  the  scattered  waves,  k  <.01  say,  then  the 
maximum  power  received  by  the  antenna  (at  9q  =  rr/2)  is  primarily  that 
of  the  direct  wave ;  i .  e .  , 


>  =  P  +  P  P 

T  D  s  D 

max 


9  =  j 
o  2 


A.  _2  .  2  .  nd. 

t  EDs,n  <-> 


(Eq.  A-7) 


where  P^  =  total  received  power. 


The  minimum  power  received  by  the  antenna  (at  9  =  0)  is  a  random 

o 

variable  whose  average  value  is  given  by  (A-6).  So  the  average  value 
of  the  ratio  of  maximum  power  to  minimum  power,  or  the  null  depth, 
is  given  by 


max 


min 


> 


<P  > 
s 


<e2> 

s 


2  sin 


2 


(Eq.  A-8) 


The  expression  (sin  x / 2) / ( 1  -  JQ(x)  )  is  not  a  very  sensitive  function 

of  x  and  sal  for  x  <  tt  (if  we  take  the  first  two  terms  of  the  expansion 

for  J  ). 
o 


Hence  we  can  write 


max  L 

<P - ~2 

T  .  k 

min 


r  d  1 
for  -  <j 


(Eq.  A-9) 
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That  this  is  a  reasonable  result  can  be  seen  by  noticing  that 

sin^  (  sin  9  )  «  (  ~r~  sin^  0  for  nd/X  small.  The  average  value 
X  n  X  n 

2 

of  sin  9  over  a  period  is  1/2.  So  a  uniform  angular  distribution  of 
n 

scatterers  would  give  at  the  output  of  the  two-element  antenna  about 
one -half  the  power  output  from  an  isotropic  antenna  of  the  same  sen¬ 
sitivity. 

Expression  (A-9)  gives  the  null  depth  of  the  two-element  antenna  as  a 
2  . 

function  of  k  .  The  accuracy  of  D/F  measurement  can  be  simply 
related  to  the  null  depth  by  examining  the  pattern  of  the  direction- 

2 

finder  (e.g.  ,  Figure  A-l).  Therefore,  it  is  reasonable  to  suggest  k 

as  a  one-parameter  criterion  to  employ  in  judging  the  suitability  of  a 

site  for  direction  finding.  Its  usefulness  may  be  appreciated  by  noting 

that  Equation  A-9  is  independent  of  the  aperture  of  the  antenna,  for 

apertures  less  than  one -half  wavelength.  This  implies  that  the  choice 

of  aperture  for  narrow  aperture  antennas  is  not  determined  by  the 

accuracy  requirements  of  the  system.  Problems  of  sensitivity  and 

antenna  balance  are  then  the  primary  considerations  in  the  aperture 

2 

size  chosen.  The  parameter  k  essentially  sets  a  limit  to  the  accu¬ 
racy  obtainable  from  a  narrow  aperture  antenna  on  a  given  site. 
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Appendix  B 

The  statistical  test  employed  to  verify  that  the  observed  amplitude  has 

the  hypothesized  distribution  specifies  that  statistically  independent 

observations  be  used.  Hence,  both  the  computer  data  and  the  measured 

data  must  be  processed  to  provide  independent  points.  It  will  be  re- 

2 

called  that  for  small  values  of  k  the  Rice  distribution  is  very  closely 

approximated  by  a  Gaussian  distribution.  Initial  calculations  showed 

2 

that  the  computed  and  measured  values  of  k  would  be  less  than  .  05, 

so  if  the  data  are  distributed  according  to  the  hypothesized  density 

function,  then  they  are  almost  Gaussian  distributed,  with  a  mean  of  1 

2 

and  a  variance  of  k  /2.  Now  let  r&  be  a  sample  data  point  and  r^, 

another  sample  a  fixed  interval  (in  angle  or  frequency)  distant  from  r  . 

a 

Further,  let  r^  and  r^  range  over  the  total  number  of  sample  points. 

(The  number  of  possible  pairs  (r  ,  r. )  is  the  total  number  of  points 

a  b 

divided  by  the  interval  between  them. )  Then  the  correlation  coefficient 

of  r  and  r,  is 
a  b 


D,b  =  <(ra  -  l)(rb  '  1):> 


(Eq.  B-l) 


where  the  brackets  (<  >)  denote  an  average  taken  over  all  possible 

sample  points  (in  angle  or  frequency).  When  p  ,  =  0,  r  and  r,  are 

a  b  a  b 

uncorrelated.  Assuming  r  and  r  are  each  Gaussian  random  vari- 

21  D 

ables,  then  r  and  r  are  statistically  independent  for  intervals  where 

p&k  =  0.  Statistically  independent  samples  of  the  measured  quantity 

r  will  be  given  by  r  ,  i  =  1,  .  .  . ,  N-ri,  where  N  is  the  total  number 

i+n 

of  samples  and  n  is  the  smallest  number  for  which 


<(r.  -  l)(r.  -  1)>  «  0 

i  i+n 
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In  practice,  the  computer  program  processing  the  data  computes  p 

ab 

for  n  =  0,  1, 2,  .  .  .  and  stops  where  p  changes  sign.  This  interval 

ab 

is  then  used  to  choose  the  proper  sample  points  for  analysis. 


Determining  whether  the  sample  points  obtained  are  Rice  distributed 

2 

is  the  next  problem.  The  technique  to  be  employed  is  the  modified  X 

1  2 

minimum  method  described  by  Cramer  and  applied  by  Siddiqui  to 
radio  propagation  problems.  This  method  involves  setting  up  a  test  of 
the  hypothesis  that  the  data  is  Rice  distributed.  The  test  essentially 
states  that  if  the  hypothesis  is  true,  then  it  is  virtually  impossible  that 
a  sample  set  of  data  points  would  be  distributed  very  differently  from 
a  Rice  distribution.  The  method  proceeds  as  follows: 

2 

1.  A  maximum  likelihood  estimator  for  k  must  be  found.  Let 

r  ,  r . r_  be  M  independent  observations  of  the  variable 

1  c  M 

r.  Since  the  hypothesized  distribution  is  approximated,  for 

2 

small  values  of  k  and  r  near  unity,  by 


(Eq.  B-2) 


the  likelihood  function  associated  with  the  density  distribution 
p(r)  is 


ilO 


A  maximum  likelihood  estimator  for  k2  is  then  given  by  a 
non-constant  (i.  e.  ,  dependent  on  . rM)  solution  of 

d(k2}  l0g  L(ri . rM:k2)  =  0  (Eq.  B-4) 

From  (B-3) 


d  log  L 

d(k2) 


I''*- 


1)  /  (Eq.  B-5) 


From  (B-5)  therefore 


k2 


=  M  l  ,Fi  -  ^ 


(Eq.  B-6) 


which  is  twice  the  variance  of  the  observations.  This  is  what 
one  would  expect,  since  the  approximating  distribution  (B-2) 

is  Gaussian  with  a  variance  kZ/2. 

The  distribution  of  the  measured  data  is  compared  with  the 
distribution  predicted  by  (B-2)  with  k2  given  by  (B-6).  An 
integer  m,  is  chosen,  not  too  small  (at  least  5)  such  that  M/m 
is  greater  than  5,  preferably  greater  than  10.  Numbers  x 
X2’  '  ”  ’  Xm  -  1  are  deterinined  such  that 


Prlr<,).^p|r)dr=i  i.-, . m  -  1 

(Eq.  3-7) 


The  range  (0,®)  of  r  is  then  divided  into  m  nonoverlapping 

intervals  I  =  (0,x  ),I  =  (x,x  ),...,  I  =  (x  . ,»),  such 

1  12  12  m  m  -  1 

that  the  expected  number  of  observations  in  each  interval  is 

M/m.  Let  the  actual  number  of  observations  falling  in  these 

intervals  be  f . f  .  Then 

1  m 


m( 

f  \2 

f  -  M) 

2  r->  \ 

,  i  ml 

■  -  — j 

I «? 

\  i=l  / 

o 

M 

m 

M| 

(Eq.  B-8) 


is  approximately  a  y  variate  with  m-2  degrees  of  freedom. 

A  critical  probability  level  a  (=  .  05  or  .  01,  say)  is  assigned 
2 

and  v  is  defined  to  be  the  number  such  that 

Pr(x2  *  X2)  =  a  (£q-  B-9) 

2  2 

If  the  observed  value  yQ  <  y^  the  hypothesis  that  (B-2)  is 

the  density  distribution  function  for  r  is  accepted.  If 
2  2 

y ^  >  y^  the  hypothesis  is  rejected.  In  the  present  work, 
the  number,  M,  of  independent  observations  ranged  between 
35  and  70,  so  m  was  chosen  to  be  equal  to  6.  The  calcula¬ 
tion  of  the  x.  explicitly  from  (B-7)  is  quite  difficult,  but 
3  1 

Norton  has  derived  an  approximation  for  small  values  of  k. 

Since  r  is  approximately  normally  distributed  for  small 

values  of  k,  it  is  appropriate  to  attempt  a  Taylor  series 

expansion  of  the  integral  (B-7)  in  powers  of  k.  Such  an 

expansion  provides  an  explicit  expression  for  the  quantities 

x.  associated  with  P  (r  <  x.)  =  i/m. 
l  r  i 


I 

I 

I 

I 

I 

I 

l 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 


112 


n 


! 


i 


X.  sa 

1 


1  + 
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(Eq.  B  - 10) 


Here  y  is  the  value  of  the  argument  of  the  normal  probability 
function  (i.  e.  ,  the  limiting  case  as  k  -*  0)  for  the  case 


P  (r  <y  .) 
r  oi 


2 

exp  (-  7  )  dx  =  — 
c  m 


(Eq.  B-ll) 


For  example,  y  =  -  .9672  for  P  (r  <y  .)  =  1/6  =  .  1667; 
ol  r  ol 

y  =  0  for  P  (r  <  y  )  =  3/6  =  .  5000;  y  =  .  4350  for  P  (r  <  y  J 
oi  r  o3  o4  r  o4 

=  4/6  =  . 6667. 


An  example  of  this  procedure  follows.  Thirty-six  statistically  inde¬ 
pendent  samples  of  r  are 


.993 

.943 

.  946 

.  944 

.  963 

.  972 

1.017 

.  958 

1.038 

1.043 

1. 039 

.995 

.989 

1.019 

.974 

.999 

1. 005 

1. 021 

1.022 

1. 029 

1.016 

1.003 

1. 061 

.967 

1.043 

1. 025 

.  980 

1.002 

.983 

.948 

.  .979 

1. 012 

.  964 

1. 001 

1. 007 

1. 055 

From  these  samples 

k2  =  .0021 


r 

i ! 
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The  resulting  class  intervals  and  the  observed  frequencies  are 


Class  Intervals 

0  -  .969 
.  969  -  . 986 
. 986  -  1. 000 
1. 000  -  1. 014 
1. 014  -  1. 031 
1. 031  -  * 


F  requencies 

8 

5 
4 

6 
7 
6 


The  expected  frequency  in  each  interval  is  6.  So  we  have 


1. 667 


If  we  choose  a  =  .  05,  for  four  degrees  of  freedom  the  critical  value  of 
2  . 

Xa  1S  9-488.  Hence,  the  data  are  consistent  with  the  hypothesized 

distribution.  In  fact,  the  probability  that  such  a  sample,  as  measured 
2 

by  y  ,  comes  from  the  hypothesized  distribution  is  more  than  70  per 
cent. 
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Appendix  C 

2 

The  fact  that  a  single  measurement  of  k  (frequency)  at  one  angle  is 
a  good  (within  a  factor  of  two)  estimate  of  the  average  value  of  these 
quantities  appears  to  have  some  relationship  to  the  familiar  concept 
of  ergodicity.  Briefly,  an  ensemble  of  time  series  is  ergodic  when 
two  conditions  hold:  1)  the  ensemble  is  stationary  in  the  strict  sense; 

2)  no  strictly  stationary  subset  of  the  ensemble  has  measure  other 
than  0  or  1.  For  ergodic  ensembles,  the  following  theorem  is  true; 
the  average  value  of  any  function  of  the  random  variable  is  the  same 
whether  the  average  is  taken  over  the  entire  ensemble  or  over  time 
for  a  particular  member  of  the  ensemble.  Typically,  in  communica¬ 
tion  problems,  very  few,  usually  only  one,  members  of  the  ensemble 
are  available  for  analysis,  so  the  "ergodic  hypothesis"  is  made  that 
the  time  series  being  analyzed  is  a  member  of  an  ergodic  ensemble. 
Thus  the  time  series  available  is  considered  to  be  representative  of 
the  entire  ensemble  and  its  characteristics  are  then  the  same  as  those 
of  the  other  members. 

In  the  calculations  made  in  Chapter  3B,  all  the  numbers  of  the  ensem¬ 
ble  are  available.  In  other  words,  the  variation  of  amplitude  with 
frequency  can  be  calculated  for  as  many  different  angles -of-arrival  as 
desired.  Hence,  the  probability  distribution  of  amplitude  across  the 
ensemble  can  be  determined  for  various  values  of  frequency.  By  anal¬ 
ogy  with  the  time-series  case,  if  it  is  the  same  for  all  frequencies, 
the  ensemble  is  stationary  in  the  strict  sense  and  the  "ergodic  theorem" 
states  that  a  single  member  of  the  ensemble  is  typical  of  all  other 
members.  The  average  value  of  any  function  of  the  amplitude  may  be 
determined  by  averaging  over  frequency  for  a  single  angle-of-arrival. 
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The  fact  that  k  (frequency)  for  a  single  angle -of-arrival  is  a  very  good 

2 

estimate  of  the  average  k  (frequency)  over  all  angles  of  arrival  indi¬ 
cates  that  the  conclusion  of  the  ergodic  theorem  holds  for  these  en¬ 
sembles.  The  conditions  of  the  theorem,  however,  are  not  true  for 

all.  For  example,  in  the  two  configurations  examined  most  closely 

2 

(Tables  3-1  and  3-2),  the  distribution  using  the  value  of  k  (angle) 

averaged  over  the  ensemble  for  all  frequencies  described  the  data 

2 

(within  the  5  per  cent  criterion  for  (r  probability)  for  more  than  95 
per  cent  of  the  separate  frequencies  for  which  the  calculations  were 
made.  The  fact  that  this  distribution  did  not  hold  100  per  cent  of  the 
time  indicates  that  it  is  highly  probable,  but  by  no  means  certain, 
that  the  ergodic  theorem  can  be  used  as  a  justification  for  using  a 
single  member  of  the  ensemble  as  representative  of  all. 
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Appendix  D 


The  material  presented  in  this  Appendix  is  a  previous  alternate  ap¬ 
proach  to  direction -finding  in  the  presence  of  reradiators  at  the  site. 
The  method  suggested  is  a  means  of  resolving  two  or  more  compo¬ 
nents  of  a  multipath  field  with  a  broad -beamed  antenna.  The  method, 
as  proposed,  has  two  significant  weaknesses: 

1.  In  order  to  resolve  several  adjacent  sources  highly 
accurate  knowledge  is  required  of  the  single-source 
antenna  pattern  and  the  measured  multiple -source 
antenna  pattern. 

2.  All  reflectors  are  assumed  small  with  respect  to  a 
wavelength;  i.  e.  ,  they  are  assumed  to  be  point 
sources  of  reradiation. 

The  second  weakness  in  this  approach  led  to  the  statistical  procedure 
described  in  the  body  of  this  work. 

The  problem  under  consideration  is  to  determine  the  angular  position 
of  each  of  several  sources  (i,  e.  ,  the  target  transmitter  and  various 
reflectors,  the  transmitter  position  is  fixed)  when  all  are  radiating 
sine  waves  of  the  same  frequency  but  of  different  amplitudes  and 
phases.  The  configuration  to  be  investigated  is  shown  in  Figure  D-l. 


1  17 


0 


a2**2 


a3**3 


FIGURE  D-l  MULTIPLE  SOURCE  DISCRIMINATION 


a, a, a,  ...  are  the  amplitudes  of  the  sources ,  and  are  the 

l  Lt  J  1  u 

phases  of  the  sources,  upon  arrival  at  point  P.  Thus  .  .  .  are 

affected  not  only  by  the  relative  phase  differences  between  each  of  the 

sources,  but  also  by  the  differences  in  path  lengths  fro  n  each  source 

4 

to  P.  The  phase  angles  are  all  given  with  respect  to  a  phase  reference 
at  P. 

Two  methods  of  determining  0  ,  9  ,  9^,  .  .  .  suggeft  themselves.  One 
is  to  move  the  point  P  and  investigate  a  segment  of  the  radiation 
field  of  the  source  distribution.  If  a  large  enough  sample  of  the  radia¬ 
tion  field  is  chosen,  it  should  be  possible  to  calcualte  what  source 
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distribution  produced  it.  This  method,  however,  presents  practical 
difficulties  in  that  an  adequate  sample  may  be  many  wavelengths  long 
and  at  low  frequencies  may  require  motion  over  considerable  dis¬ 
tances. 

The  other  method  will  be  explored  here.  It  involves  placing  a  direc¬ 
tional  antenna  at  P,  rotating  it  about  its  phase  center,  and  examining 
the  output  as  a  function  of  aiming  angle.  From  knowledge  of  the  an¬ 
tenna's  response  to  a  point  source  and  of  the  measured  pattern,  it  is 
possible  to  deduce  what  the  source  distribution  is.  This  may  be  done 
in  either  of  two  ways: 

1.  The  response  of  the  antenna  to  a  great  variety  of  source 
distributions,  amplitudes,  and  phase  angles  may  be  cal¬ 
culated  and  the  measured  response  compared  to  each  of 
these.  The  calculated  response  most  closely  corresponding 
to  the  measured  response  will  give  the  source  distribution. 

2.  An  analysis  of  the  measured  pattern  based  on  its  Fourier 
expansion  will  give  directly  the  values  of  0^ ,  0^,  .  .  .  This 
analysis  is  presented  in  this  Appendix. 

It  should  be  emphasized  that  the  analysis  which  follows  is  superfluous 
if  the  antenna  placed  at  P  has  a  sufficiently  large  aperture.  For  then 
the  antenna  will  possess  sufficient  resolution  to  distinguish  between 
various  sources.  The  following  calculations  are  necessary  for  large 
beamwidth  antennas  only. 

Now  assume  the  antenna  at  P  has  a  field  response  of  A(0).  Then  the 
response  of  the  antenna  to  a  unit  magnitude  point  source  at  @  is 


A(0  -  9j).  Hence,  the  signal  received  by  the  antenna  placed  at  P  of 
Figure  D-l  can  be  represented  in  the  phasor  diagram  shown  in  Figure 


3 


D-2. 


FIGURE  D-2  ADDITION  IN  THE  COMPLEX  PLANE 


The  orientation  of  the  x  and  y  axes  is  chosen  arbitrarily  as  a  phase 
reference.  The  x-component  of  the  resultant  phasor  is 

n 

R  (9)  =  R(9)  cos  n  (9  )  =  )  a.  A(9  -  9.)  cos  $.  (Eq.  D-l) 

x  L,  iii 

i  =  1 

and  the  y-component  of  the  resultant  phasor  is 

n 

R  (9)  =  R(9)  sin  n(9)  =  >  a.  A(9  -  9.)  sin  $.  (Eq.  D-2) 

y  L>  iii 

i  =  1 
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In  complex  notation,  Equations  D-l  and  D-2  become 


n  n 

J  * 

R  (9  )  =  R<9)  ejp(9)  =  V  a.  A(9  -9.)  e  1  =  V  a.  A(9  -  9  ) 

L  i  i  L>  i  i 

i  =  1  i=  1 

(Eq.  D-3) 

Note  that  R(0)  is  a  measurable  quantity  since  R(9)  and  p(9)  are  meas¬ 
urable. 


Clearly  R(9)  is  a  periodic  quantity,  with  period  =  2rr.  Hence  it  may  be 
expanded  in  a  Fourier  series: 


03 

R(9>  =  l  ?kJk0 

k  =  -os 

where 


(Eq.  D-4) 


R(9  )  e'^9  d8 


(Eq.  D-5) 


We  know  that  A(9)  is  periodic  also,  hence  it  has  the  Fourier 
expansion 


A(9) 


(Eq.  D-6) 
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where 


°k  =  2tt 


C 

ij 


A(9  )  e"'ik9d0 


(Eq.  D-7) 


Hence,  the  right-hand  side  of  Equation  D-3  may  be  written 


"  "  „  jk(9  -  9.) 

I  *1  A(9-9i» =  I  »i  I  %e 


i  =  1  k  =  -“ 


.  y  o\  ejk5  V  a. 
L  k  L  i 


-  jk  9 . 


(Eq.  D-8) 


k=-°° 


Finally,  from  Equations  D-4  and  D-8,  we  have 


0,  A  -ikq- 

k  ~  V  ~  1 

£7=  6k  =  L  aie 

k  ;=i 


(Eq.  D-9) 


Equation  D-9  represents  a  set  of  nonlinear  simultaneous  equations,  one 

for  each  value  of  k.  The  left-hand  side  of  Equation  D-9  is  known,  since 

it  is  computed  from  the  known  quantities  A(9)  and  R(8).  Since  there 

are  2n  unknown  quantities,  a,  ...  a  and  0  ....  9  ,  we  expect  that 

lnln 

only  2n  of  Equation  D-9  will  be  independent,  and  we  can  choose  the 
most  convenient  set  for  the  solution  of  the  problem. 

If  we  pick  k  =  +  (n-1)  ,  ...  +1,  0,  -1,  ...,  -n,  and  make  the  substi  - 
.  o  +  i  0 ; 

tutions  S.  =  e  we  get  the  equations 


. 


~  ~  -(n-1)  ~  c  "(n-1)  ,  ~  c  -(n-1) 

o  ,  =  a  S  +  a  b  ,  +  a  S 

n-1  11  2  2  mn 


~  ~  -  -(n-2)  ,  ~  Q  "(n_2)  ,  ~  q  -(n"2) 

5  =  a  S  +  a.S  .  +  a  S 

n-2  11  22  m  n 


,  .  .  .  ,  +  a 


~  ~  n-1  ~  c  n-1  ,  ~  c  n-1 

-(n-1)  11  22  nn 


-  n 
=  a  S 
i  i 


~  n 
+  a_S0 
2  2 


.  +  a  Sn  (Eq.  D - 1 0) 

n  n  n 


This  set  of  equations  can  be  solved  by  the  following  procedure: 

a.  Multiply  the  equation  for  k  =  +(n-l)  by  and  subtract  it  from 
the  equation  for  k  =  +(n-2)  to  form  the  first  equation  of  a  new 


b.  Multiply  the  equation  for  k  =  +(n-2)  by  and  subtract  it  from 
the  equation  for  k  =  +(n-3)  to  obtain  the  second  equation  of  a 
new  set. 

c.  Repeat  the  process  for  (2n-l)  pairs  of  equations.  There  will 
result  (2n-l)  new  equations  with  a^  eliminated. 

d.  Repeat  the  entire  process  on  the  new  set  of  (2n-l)  equations, 
multiplying  by  to  obtain  a  third  set  of  (2n-2)  equations 
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with  a^  and  a^  eliminated.  The  process  is  continued  until  all 

the  a.'s  are  eliminated. 

1 

The  procedure  results  in  the  set  of  n  equations 


T  +  b  ,  nu,  +  6  ,  ,.U, . +  6  U  =  0 

-n  - (n - 1 )  1  -(n-2)  2  on 


T  ,  n  +  6  ,.u,  +  6  ,  ,.u_,  . . . ,  +  7.  ,u  =  o 

-(n-1)  -(n-2)  1  -(n  -  3)  2  +1  n 


6 


-1 


where 


■  6  U,  +  6  ,  ,U, 

.  .  .  ,  +  s'  ,  l4U  =0 

o  1  +12 

+(n- 1)  n 

(Eq.  D 

u,  .  -(s1+s2  +  s3  +  s4. 

....  +sn) 

U  =  +(S  s  +  s  s . 

2  'i2  23 

+  S  S  ) 

n-1  n 

U  =  (-l)n  (S  S  ...  S  ) 
n  \  c  n 


(Eq.  D-  12) 


From  Equations  D-ll  it  is  possible  to  solve  directly  for  the  IF.  Once 
the  U.  are  known,  the  S.  can  be  computed  by  the  following  technique: 

a.  Add  Sj  to  both  side  of  the  first  of  Equations  D-12. 

b.  Multiply  each  of  Equations  D-12  by  m,  where  m  is  the  nu 
merical  subscript  of  U. 
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c.  Add  the  resulting  equations  to  yield  the  following  polynomial 
S^  +  YjS?'1 . +  Yn^Sj  +  Yn  =  0  (Eq.  D-13) 

<N/ 

where  y.  are  the  (complex)  solutions  for  U  arrived  at  from  Equa- 
1  1 

tions  D-  1 1 . 

The  elimination  procedure  leading  to  Equation  D-3  might  just  as  easily 
have  used  any  of  the  S. .  From  this,  it  follows  that  the  n  complex 
roots  of  Equation  D- 3  are  the  values  of  S.,  i  =  l,  n.  Having  found 

S  ,  S  ,  ....  S  (the  numbering  is  obviously  arbitrary)  we  can  substi- 

X  L*  IX 

tute  into  any  n  of  Equations  D-10  to  find  the  a,,  a_ . a  associ- 

1  L  n 

ated  with  each  of  them.  The  amplitude,  phase,  and  angular  position  of 
each  source  of  Figure  D-l  have  now  been  completely  determined. 

The  choice  of  +(n-l) . +1,  0,  -1,  .  .  .  ,  -n  for  the  values  of  k  pre¬ 

ceding  Equations  D-10  is  the  result  of  the  desire  to  require  as  few 
harmonics  as  possible  in  the  Fourier  expansion  of  the  antenna  pattern. 
Thus  it  is  necessary  to  have  at  least  as  many  harmonics  as  sources. 

An  interesting  consequence  of  this  statement  is  that  the  ability  of  an  an¬ 
tenna  satisfying  this  criterion  to  distinguish  between  various  sources 
is  limited  only  by  the  accuracy  with  which  R(0)  and  A(8)  can  be  meas¬ 
ured.  The  solution  presented  above  is  the  same  regardless  of  whether 
two  sources  are  lu  apart  or  90°  apart  in  angular  position.  Thus,  if  the 
initial  measurement  were  exact,  the  antenna  system  described  (antenna 
plus  computation  equipment)  would,  in  effect,  possess  infinite  resolving 


power  . 


As  an  example  of  this  procedure,  suppose  that  the  antenna  pattern  is 

A(9  )  =  1  +  cos  9  +  cos  29  (Eq.  D-14) 

The  half-power  beamwidth  is  72°.  By  the  method  presented  in  this 
appendix,  such  an  antenna  should  be  capable  of  resolving  two  sources 

/■v 

a  and  a.  at  angular  positions  9^  and  9  respectively.  In  other 
words,  Equation  D-13  should  reduce  to  an  identity  for  S  =  eJ  1  or 
e^2.  From  D-14  we  have 


a„  =  a. 


\_ 

2 


a ,  =  a 


1 

2 


(Eq.  D-15) 


a  -  1 

o 


The  response  of  the  antenna  to  the  two  sources  is 


R(9)  =  1  +  cij  cos  (9  -  9j)  +  a2  cos  (9  -  9  2) 


+  a.j  cos  2  (9  -  9j)  +  a2  cos  2  (9  -9^) 


F rom  which 


y  1  /-  2j9l  .  ~  2j92\ 

6-2  =  llV  +  V  ) 

y  1  f~  j°l  ~  2j92l 

S-l  =  l(V  +  V  ) 


3  =  a  +  a 

o  1  2 


(Eq.  D-16) 
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From  D-15  we  have 


(Eq.  D-17) 


(Eq.  D-18) 


which,  following  D-10,  are  all  we  need. 


The  solutions  to  D-ll  are 


L  T  2 

-  1  1  0 


(Eq.  D-19) 
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which,  when  the  values  of  6.  are  substituted,  become 
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